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This thesis presents theoretical work in s- and p-wave resonantly paired Fermi 
gases at zero temperature. In the BEG regime of the wide-resonance s-wave BGS-BEG 
crossover, the chemical potential, speed of sound, condensate depletion, and ground 
state energy are computed. The quantities are calculated in a low density approxima- 
tion. The bosonic scattering length is computed diagrammatically to be approximately 
0.60a with a the fermionic scattering length, confirming a well-known result. The per- 
turbative approach breaks down in the intermediate BGS-BEG crossover regime which 
is not investigated. Quantum corrections are computed in the low density expansion, 
and the mass-imbalanced two-species Fermi gas is investigated. For identical fermions, 
it is discussed how p-wave Feshbach resonances naturally fall into two classes of "weak" 
and "strong" , depending on the short-distance physics. It is shown how bound fermionic 
trimers appear in the strongly-resonant superfluid. The appearance of the trimer under 
the strengthening of the resonance is investigated, and the lifetime of p-wave diatomic 
molecules due to coUisional relaxation into trimers is estimated. 
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Chapter 1 



Introduction 



The history of superfluidity dates back to the discovery of superconductivity 
(charged superfluidity) in mercury by KamerUngh Onnes [54]. Onnes observed that 
the electrical resistance of mercury disappears completely as the system is cooled below 
a critical temperature. A detailed microscopic explanation of the phenomenon of super- 
conductivity proved elusive until the breakthrough in the 1950's by Bardeen, Cooper, 
and Schrieffer [5, 6]. The work of these authors (BCS theory) described how an arbi- 
trarily weak attractive interaction in a degenerate Fermi gas leads to the formation of 
Cooper pairs [18], composed of weakly correlated pairs of fermions. The BCS theory 
explained the resulting appearance of a minimum excitation energy, or energy gap, and 
proved hugely succesful in understanding properties of superconductors. In particular, 
the BCS theory was able to correctly describe superfluid ^He (fermionic atoms) but not 
superfluidity in liquid "^He (bosonic atoms). 

The development of the theory of Bose-Einstein condensation (BEC) occurred 
in parallel. In 1925 Einstein predicted [26] the macroscopic occupation of the single- 
particle ground state below a critical temperature in systems of particles obeying Bose 
statistics. This followed ideas by Bose on the quantum statistics of phonons [11]. In 1938 
it was proposed by London [51, 52] that superfluidity in ^He appears as a consequence 
of Bose-Einstein condensation. However, liquid ^He is a strongly interacting system, 
very different from the weakly interacting limit studied by Einstein. In particular, the 



number of particles occupying the single-particle ground state is strongly suppressed 
even at zero temperature, making the verification of Bose-Einstein condensation in 
liquid ^He extremely difficult. 

Considerable experimental effort was spent on attempting to obtain a BEC in a 
weakly interacting system. This was first achieved in 1995 in dilute gases of alkali atoms 
consisting of rubidium [3], sodium [21], and lithium [13]. In these experiments, the gas 
was cooled to temperatures for which the thermal de Broglie wavelength was of the order 
of the interparticle spacing and the BEC of weakly interacting atoms could be observed 
very clearly from the bimodal profile of the momentum distribution of the gas below 
the critical temperature. For a review of the subject of Bose-Einstein condensation in 
dilute gases see e.g. Refs. [19, 48, 55]. 

Soon after the observation of BEC in dilute gases of bosonic atoms, considerable 
attention was shifted to the study of Fermi gases. Applying the same experimental 
methods as in the bosonic gases proved challenging, as the s-wave interactions necessary 
for evaporative cooling are absent in spin-polarized Fermi gases. Thus the cooling of 
Fermi gases requires the presence of two distinguishable atoms, either in the form of 
two hyperfine spin components of the same atom or in mixtures of two atomic species. 
The first realization of a quantum degenerate two-component Fermi gas was obtained in 
a system of ^"^K [23] where a temperature of T ~ 0.2Tp was achieved {Tp is the Fermi 
temperature). Soon experiments followed in which quantum degeneracy was observed 
in ®Li [72, 66] using sympathetic cooling with the '^Li isotope. 

In two-component ultracold Fermi gases, it should be possible to observe super- 
fluidity according to the BCS theory. However, in these first experiments which reached 
quantum degeneracy, the critical temperature required to enter the superfluid phase 
was very small due to the diluteness and weak interactions in these gases. It was soon 
recognized [37, 70] that it is possible to enhance the weak interactions by the use of a 
scattering resonance known as a Feshbach resonance [27], thereby increasing the critical 



temperature. 

The use of Feshbach resonances opened up the possibihty of studying both BCS 
superfluidity and BEC in the same system. These resonances are a characteristic of 
the two-body interactions and they aUow for the precise control of the strength of 
interparticle interactions through the tuning of an applied magnetic field. For a weak 
and attractive interaction in a two-component Fermi gas, BCS theory is valid and the 
system is expected to form a BCS superfluid below a critical temperature. If the two- 
body interaction is made progressively stronger by the tuning of a magnetic field, at 
some critical strength of the interaction a two-body bound state becomes possible. The 
bound states are then bosons and form a BEC at sufficiently low temperatures. A 
BEC in a two-component Fermi gas was first observed in 2003 through the typical 
bimodal distributions [33, 40, 78]. The observation of the BCS superfluid phase is more 
subtle as the normal and superfluid phases share many properties. Convincing proof 
of superfluidity was obtained in 2005 through the observation of quantized vortices on 
both sides of a Feshbach resonance [77] . 

The possibility of a BCS-BEC crossover was originally proposed by Eagles in 
1969 [24], and later by Leggett [47], and by Nozieres and Schmitt-Rink [53]. These 
authors extended the mean field BCS model, valid in the weakly interacting BCS regime, 
to an arbitrary interaction strength, i.e. to well outside the regime of the original 
model's validity. It was discovered that under the strengthening of the interaction, the 
BCS gap equation, the mean field theory which describes the BCS state, evolves into 
the Schrodinger equation of a bound pair of fermionic atoms, correctly predicting the 
binding energy of the diatomic molecule. Furthermore, the results of these authors 
suggested that no phase transition occurs in the crossover, and thus the BCS and BEC 
phases are qualitatively the same phase. These results connected the two types of paired 
superfluidity, within mean field theory, and demonstrated how BCS and BEC are two 
extreme limits of the same system. 



Despite the success of the gap equation in describing the BCS superfluid and the 
bound pairs in the BEC regime, it remained an open question whether it could correctly 
predict properties of the many-body system on the BEC side of a Feshbach resonance. 
The short answer is that it cannot; the gap equation, when applied to the BEC side 
of the crossover, implicitly assumes that the perturbative parameter is the interaction 
strength between the Feshbach molecules, i.e. that interactions are weak and scattering 
proceeds in the Born approximation. However, the Born approximation is known to 
fail to describe interactions between molecules properly [59]. It follows that while the 
BCS-BEC gap equation may describe the physics of the BEC regime at a qualitative 
level due to the lack of a phase transition, it is not to be trusted quantitatively. Thus 
there is a need for a reliable technique in the BEC regime, beyond mean field theory. 
This is one of the main subjects of this thesis (see also Ref. [50]). 

In this thesis, the small parameter used in the BEC regime is the gas parameter 
na^. Here, n is the density of fermions and a the s-wave scattering length between 
the two components of the Fermi gas. In the gas parameter expansion the chemical 
potential, ground state energy, speed of sound, and condensate depletion are computed. 
The results are seen to match the results of the standard dilute Bose gas [1, 29], how- 
ever with the boson-boson scattering length related to the underlying fermion-fermion 
scattering length by 

Ob « 0.60a. (1.1) 

This result was first computed in Ref. [59] using a coordinate space formalism and 
the result is confirmed in this thesis using a diagrammatic technique. While the work 
reported in Ref. [50] and also in this thesis was in progress, the result (1.1) was computed 
independently by Brodsky et al [14] using essentially the same diagrammatic technique. 
At the crossover, the scattering length a diverges. Thus, the parameter na^ is 
no longer small, and properties of the gas close to the Feshbach resonance will not be 



computed. 

In the 1950's, perturbative corrections in powers of no? to the chemical potential 
and ground state energy were computed in the standard dilute Bose gas [45, 46, 76, 65, 
38] . It is explicitly demonstrated in this thesis how the fact that the bosons are composed 
of fermions does not change the first two such corrections, apart from replacing Cb by 
0.60a. Even higher order contributions will depend explicitly on the presence of fermions 
in the problem but will not be considered here. 

An interesting system is a two-component Fermi gas where the two components 
have different masses. For small mass imbalances, this system still displays the usual 
BCS-BEC crossover physics, with the few-body couplings slightly changed. However, 
as the mass imbalance is increased, the system becomes less stable until finally above 
a mass-ratio of 13.6 between the fermion components, the three-body system has been 
found [57] to contain bound states similar to the Efimov states [25] in systems of bosons. 
Properties of the mass imbalanced system such as the atom-molecule and molecule- 
molecule scattering lengths and the loss rate due to coUisional relaxation have been 
calculated in Refs. [56, 57, 59] by use of a real-space technique, and the results will be 
confirmed in this thesis using a diagrammatic technique. 

To date, experiments in ultracold gases of fermionic atoms have primarily studied 
gases consisting of two hyperfine states of the same atom. However, recently there has 
been a considerable interest in gases of identical fermionic atoms interacting close to a 
p-wave Feshbach resonance. The observation of p-wave Feshbach resonances has been 
reported in ^"^K [69] and ^Li [67], and theoretical work on the subject includes Refs. 
[36, 22, 35, 17, 34, 49]. In the p-wave superfluid, the BCS and BEC phases are no 
longer extreme limits of the same phase, rather they are separated by a genuine phase 
transition [74, 22], or in some cases a sequence of phase transitions [35, 17, 34]. The 
phase transition may be topological in some cases [62, 35, 43]. 

j»-wave Feshbach resonances naturally fall into two categories of weak and strong. 
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Mean field theory in the crossover, as developed in Refs. [35, 17, 34], is only strictly 
valid for the weak resonances. On the other hand, the p-wave resonance studied in the 
experiment of Refs. [69, 30] is strong [34]. It is therefore of interest to study the strong 
resonances. It was first noticed by Y. Castin and collaborators [16] that in this regime, 
three fermions will form a bound state with angular momentum 1. These trimers are 
very strongly bound, with a size of the order of the short range physics which led to the 
formation of the Feshbach molecule. In this thesis, the appearance of the trimer under 
the strengthening of the resonance will be studied, and the decay rate due to inelastic 
collisions in which two molecules turn into a trimer and an atom will be estimated. 

In this thesis the effects of a trapping potential will not be considered. Of course, 
the trapping potential is an essential ingredient in the experimental study of ultracold 
gases. However, the local density approximation has proven quite succesful in describing 
most properties of the system as long as the trapping potential is sufficiently smooth. 
For a review of trapped gases see e.g. Ref. [19] . 

The thesis is organized as follows. Chapter 2 presents some simple properties of 
the BCS-BEC crossover, and a model of the Feshbach resonance physics is introduced. 
Chapter 3 contains calculations of properties of the BEC-BCS crossover in the BEC 
regime, and in chapter 4 higher order corrections to the chemical potential and ground 
state energy are computed in the gas parameter expansion. The two-component mass- 
imbalanced Fermi gas in the BEC regime is studied in chapter 5, with emphasis on few 
body scattering lengths and collisional losses. The focus then shifts to p-wave resonantly 
coupled superfluids. In chapter 6 the system is described and the model introduced. 
Chapter 7 investigates the behavior of strongly-resonant p-wave superfluids and the 
appearance of bound trimer states under the strengthening of the resonance. 

For simplicity of notation, Planck's constant h and Boltzmann's constant ks have 
been set to 1. 



Chapter 2 
Physics of the BCS-BEC crossover 



Soon after the achievement of Bose-Einstein condensation in dilute gases of bosonic 
alkah atoms, considerable attention was instead focused on gases comprised by fermionic 
atoms. The difference in quantum statistics makes these systems quite distinct at ul- 
tracold temperatures. As temperature is lowered, a bosonic system goes through a 
phase-transition to a BEC and this happens even in the absence of interactions. The 
Pauli principle, on the other hand, prevents two identical fermions from occupying the 
same single-particle state. Thus, in a non-interacting Fermi gas at zero temperature, 
the fermions will instead occupy the lowest energy states with exactly one fermion in 
each state. Hence, superfluidity in a gas of fermions can only stem from the presence 
of interactions. As the temperature of a Fermi gas is decreased, the crossover from a 
classical to a quantum behavior is smooth. 

Consider a gas of spin 1/2 fermions at low temperature. If the interactions in the 
spin singlet channel are arbitrarily weak and attractive, the Fermi sea will be destabi- 
lized towards the formation of Cooper pairs [18]. The pairing occurs between fermions 
with equal and opposite momenta, the pairs being weakly correlated in momentum 
space. The pairs are highly overlapping in real space and should not be thought of 
as molecules. The ground state of the system is then no longer a Fermi sea, rather 
it is a superposition of Cooper pairs, the primary contribution arising from fermions 
close to the Fermi surface. This is explained in the famous work of Bardeen, Cooper 





Figure 2.1: An illustration of the BCS-BEC crossover. The BCS regime on the left is 
illustrated as the pairing in momentum space of fermions, while the EEC regime on the 
right is illustrated as the real-space pairing of atoms into molecules. Paired fermions 
in the central regime between the two limits share features with the paired fermions in 
both limits. 



and Schrieffer [5, 6]. The system studied by these authors consisted of pairs of spin- up 
and spin-down electrons and the resulting many-body system was a superconductor (a 
charged superfluid). In the context of ultracold two-species Fermi gases, the regime 
described is known as the BCS regime. 

Imagine then that by some means the attractive interaction is increased. As the 
strength of the interaction is increased, a bound molecular state between two fermions 
becomes possible. The resulting bound state is then a boson, and a gas of these diatomic 
molecules may form a EEC. This is the essence of the BCS-BEC crossover; as the 
interaction between fermions is increased, there will be a continuous crossover from a 
BCS superfluid to a BEC of diatomic molecules. Thus the BCS and BEC regimes are 
merely limiting cases of the same system. The crossover is illustrated in Fig. 2.1. In the 
central region between the two limiting regimes, the bosonic pairs retain some properties 
of both Cooper pairs and diatomic molecules. 

With the advent of BCS theory in the 1950's, the limiting regimes in the BCS- 
BEC crossover were well understood. However, the connection between these types of 
superfluids was not appreciated until the work of Eagles in 1969 [24], and the later 
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works of Leggett [47] as well as Nozieres and Schmitt-Rink [53] in the 1980's. These 
authors started from a BCS mean-field theory which is valid only for weak attractive 
interactions. They extended the original model to arbitrarily strong attractive interac- 
tions, well outside the range of validity of the model. The mean field equation which 
describes properties of the BCS superfluid is known as the gap equation and it was dis- 
covered [47, 53] that the gap equation evolves into the Schrodinger equation of a bound 
pair of fermions under the strengthening of the interaction. Thus the results of these 
authors connected the BCS superfluid with the BEC of condensed pairs of fermionic 
atoms, within a mean field approach. Finite temperature studies [53, 64] further showed 
that the transition temperature evolves smoothly between the two regimes. The results 
suggested a smooth crossover without any phase transition. It should be emphasized 
that the results obtained by this method, while qualitatively correct due to the lack of 
a phase transition, are not quantitatively correct in the BEC regime, since the mean 
field approach is not strictly valid. This point will be discussed further in chapter 3. 

The crucial tool which allows the experimental control over the interaction strength 
in a two-component Fermi gas is a Feshbach resonance [27, 28, 71]. A Feshbach res- 
onance is an intrinsic two-body phenomenon. Consider the interaction between two 
fermionic atoms; at short distances the interaction is repulsive, while at large distances 
it is dominated by the van der Waals attraction which goes as —C^/r^. In general, the 
interatomic potentials contain several vibrational bound states. A Feshbach resonance 
occurs when the energy of a bound state in one scattering channel (usually called the 
closed channel) coincides with the energy of free particles scattering in another channel 
(the open channel). The two spin channels will in general have different magnetic mo- 
ments. Typically, the closed channel is an approximate spin singlet, the open channel 
an approximate spin triplet, and the (large-distance) splitting between the energy levels 
is provided by the Zeeman effect. The Zeeman splitting may be tuned by adjusting the 
magnitude of an external magnetic field. The situation is illustrated in Fig. 2.2. The use 
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U(r) 




Figure 2.2: A Feshbach resonance occurs as the energy of the bound state in the closed 
channel crosses the energy of free particles scattering in the open channel. 



of Feshbach resonances allows a very precise tuning of the strength of the interaction, 
which may even be changed in real time. 

The presence of the (real or virtual) bound state in the closed channel can sig- 
nificantly alter scattering in the open channel. Across a Feshbach resonance, the (open 
channel) scattering length behaves as 

AS 



a{B) = an 



1 + 



B-B^ 



(2.1) 



Here, Onr is the non-resonant scattering length, the value of the scattering length a far 
off-resonance, Bq is the magnetic field position of the resonance, and /S.B is known as 
the magnetic resonance width (defined as the difference in magnetic field, Bq — B, at 
which the scattering length becomes zero). The behavior of the scattering length is 
illustrated in Fig. 2.3. 

Below, low-energy scattering properties relevant for the present work will be re- 
viewed. A model of the resonant interactions, the so-called two-channel model, will be 
considered and its relevance to the Feshbach resonance physics described above will be 
discussed. 
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Figure 2.3: The scattering length a across a Feshbach resonance. On the left side of the 
resonance a real bound state exists. 



2.1 Low energy scattering 

This section contains a short discussion of the scattering theory relevant for this 
thesis. For elastic scattering of two particles via a centrally symmetric interaction 
potential, the scattering amplitude may be written as [44] 



f{k,e) = J2i'^£ + i)p,icose)fi{k). 



(2.2) 



e=o 



Here, \k\ = \k'\ = k is the magnitude of the incoming, zizk, and outgoing, zizk' , momenta, 
while 9 is the scattering angle, with k ■ k' = k"^ cos 6. The probability of scattering 
through the solid angle $7 is related to the scattering amplitude by da/d^ = |/p. A 
centrally symmetric interaction U{r) conserves the angular momentum (. and in this 
thesis, focus will be on s- and p-wave scattering, with 



fs{k) ^ f,=o{k), 
fp{k)^h=i{k). 



(2.3) 
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The partial wave amplitudes may also be written in the form [44] 

The functions Fi{k'^) are real and may be Taylor expanded in their arguments. Assuming 

that all the expansion coefficients take values of the order of some "natural" scale, it 

is seen that for low-energy scattering (e.g. in ultracold gases) the functions Ff{k'^) are 

well approximated by the first few terms of their Taylor expansions. In particular it is 

important to note that the scattering amplitude at low energy only depends very weakly 

on the precise details of the complicated inter- atomic interaction potential U{r). 

For s-wave scattering, keeping only one term in the Taylor expansion of F(^{k'^) 

amounts to 

fs{k) = _\ ., , (2.5) 

—a ^ — ik 

where a is the s-wave scattering length. The scattering amplitude has one pole, located 

at kpoic = +ia^^ ■ This pole only corresponds to a real bound state if a > 0, otherwise 

it corresponds to a virtual bound state [44]. For a > the real bound state has the 

energy 

where m^ is the reduced mass in the two body system. The result already nicely matches 
the Feshbach resonance picture described above. However, the approximation Eq. (2.5) 
is deficient in one respect, it lacks a true resonance, i.e. a quasistationary state with a 
finite lifetime [34]. In order to correctly capture the possible existence of such states, 
the next order expansion 

F,(fc2) = -a-i + ^rofc' (2.7) 

is needed. The parameter vq is known as the effective range of the interaction potential 

U{r). 
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An important energy scale is set by the effective range 

4 



(2.8) 



mrQ 

This energy scale is related to the magnetic field width by the approximate relationship 
To « 4m/u^anj,Ai?^ [34]. The gas contains a finite atom density n, and consequently an 
additional length scale is set by the average particle spacing n~^'^. The corresponding 
energy scale is the Fermi energy^ ep given by 

ep = 1^. (2.9) 

kp is the Fermi momentum. The two energy scales allow the definition of a dimensionless 
parameter [34] 

7. = ^J^, (2.10) 

TT V EF 

which measures the width of the resonance. Resonances for which 7s ^ 1 are termed 
narrow while those with 7^ ^ 1 are called wide. The precise criterion for the resonances 
to be wide was first derived in Ref. [15]. 

The distinction between wide and narrow resonances is quite important. Since the 
effective range only varies weakly across a Feshbach resonance, for narrow resonances 7^ 
provides a small parameter in which a quantitatively correct (in the zero width limit) 
perturbation theory may be formulated across the whole crossover [34]. In contrast, no 
such simplification occurs for wide resonances; the other dimensionless parameter which 
may be constructed from the three length scales a, tq, and n^^'^ is known as the gas 
parameter, v}i^a. Since the scattering length diverges across the resonance, no small 
parameter exists close to a wide Feshbach resonance. 

Most current experiments in ultracold two-component Fermi gases are carried out 
close to wide resonances. As these are the most relevant, the focus in the parts of this 
thesis concerning s-wave Feshbach resonances (chapters 3, 4, and 5) will be on wide 
resonances. 



^ In the absence of two-body bound states this would have been the conventional Fermi energy. 
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It is precisely for the narrow resonances that the two-body problem may display 
a true resonance (a positive energy state with a finite lifetime) [34]. Since the s-wave 
Feshbach resonances studied in this thesis are wide, these quasistationary states will 
be ignored. It is, however, interesting to note that although most experiments are 
performed close to wide resonances, since 7^ ~ n~^i^ resonances which are wide in 
current experiments may become narrow if higher densities are achieved. 

2.2 The two-channel model 

A model which reflects the physics of the interaction in a two-component Fermi 
gas interacting close to a Feshbach resonance is provided by the two-channel model. 
It was first studied by Rumer in 1959 [63] and in the context of atomic Bose-Einstein 
condensates it was studied by Timmermans et al in 1999 [71]. Without interactions, 
the model describes the open-channel atoms and closed-channel molecules separately. 
Interaction between atoms and molecules is provided by an interconversion term. 

The model is conveniently defined in terms of a functional integral 



Z = I V'4)^Vij^Vil^{Di){DbVbe'^'''. (2.11) 



^1 and -i/^j are the fermionic atoms in two different hyperfine states (for simplicity 
denoted | and j) and b is the closed-channel molecule. The action S'tc consists of three 
parts, S'tc = Sq + Sob + Sfb. Of these. So is the free action of fermions 

So= E /^^-^*^^(^| + ^£ + m)^^, (2.12) 

where ^ is the chemical potential of the fermions. The free action of the closed-channel 
molecules is 

d\dt-b(^- + ——^+2^-e,y. (2.13) 

eo is called the detuning and is the parameter which is typically controlled by a magnetic 
field in the experimental setup. That is, the detuning controls the proximity to the 
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Feshbach resonance. The interconversion term between open-channel atoms and the 
closed-channel molecule is given by 

Sp^ = -g cPxdt [hi)^'4>i + b'4>i^^] . (2.14) 



The coupling g is the probability of converting a h particle to a pair of atoms. In itself, 
the interaction occurs at a point in space and is unphysical; the action needs to be 
supplemented by a cut-off, A ~ -^ where Rf. denotes the range of the (short-ranged) 
forces. 

It is important to note that the molecule h is only the true bound state in the 
case of a vanishing coupling between atoms and molecules. Otherwise, the true bound 
state is the "bare" molecule h surrounded by a cloud of atoms. Thus, while the bare 
molecule's size is of the range of the forces which led to its creation, i?e, the size of the 
physical molecule is of the order of the atom- atom scattering length a ^ Re. In fact, 
the contribution of the bare molecule to the physical molecule effectively disappears for 
infinitely wide Feshbach resonances, where the two-channel model may be replaced by 
a one-channel model. This will be discussed in more detail in chapter 3. 

The two parameters which characterize the scattering between the atoms are the 
scattering length and the effective range. These may be calculated in the model and 
are given by [4] 

" = 7 TTY' ^^0 = 2^- (2-15) 



It is apparent that the position of the resonance is shifted from the "bare" value eo to 



the physical detuning, 



Wo = eo ^^- (2.16) 



The shift again occurs because the bare molecule does not coincide with the physical 
molecule. 
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In terms of the parameters of the model, 

7. = -77^. (2.17) 

In the two-channel model, it is again seen that 7^ is independent of detuning loq and 
scattering length a. Thus, if 7^ ^ 1 away from resonance then this is the case across 
the whole crossover. On the other hand, if rg is small then the resonance will be wide. 
In the next chapter the effective range will be taken to zero in order to model the wide 
resonances. 



Chapter 3 
The BCS-BEC crossover in the BEC regime 



In the previous chapter some basic properties of the BCS-BEC crossover were 
described. Consider a two-component Fermi gas with short-ranged interactions. If the 
interactions between the two components are weak, the fermions will pair into the usual 
Cooper pairs and the system becomes a BCS superconductor at low temperatures. If the 
strength of the interaction is then increased by some means, at some point it becomes 
possible for two fermions to form a bound molecular state. As the strength is tuned 
further away from resonance, the bosonic molecules will form a regular BEC, and it is 
properties of this BEC which will be described in the present chapter. 

The s-wave scattering length of two different fermionic atoms diverges across the 
Feshbach resonance. However, the BCS-BEC gap equation which describes the weakly 
attractive Fermi two-component Fermi gas implicitly assumes that the scattering is 
weak. Thus, even though the gap equation correctly describes the physics on the BCS 
side of the resonance there is no reason why the gap equation should be trusted on the 
BEC side of a Feshbach resonance. Yet, in Refs. [47, 53] it was discovered that under 
the strengthening of the interaction, the gap equation develops into the Schrodinger 
equation for the bound pair of fermionic atoms and correctly predicts the binding energy 
of this bound pair. Studies by these authors furthermore demonstrated that the BCS 
superfluid evolves into the BEC without undergoing a phase transition. Thus, even 
though the gap equation is not quantitatively correct on the BEC side of a resonance. 
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it still yields qualitatively correct results. 

In order to systematically investigate the BEC regime of the weakly bound bosonic 
molecules, it is necessary to identify a small parameter. This is provided by the so-called 
gas parameter n}/^a. n is the density of fermions while rib = n/2 is the density of 
bosonic molecules. The size of a molecule is roughly of order a, as discussed in chapter 

— 1/3 

2, while the average intermolecular distance is n|_, . Thus, even though the molecules 
are rather large compared to the short-distance physics which led to their formation, 
for small values of the gas parameter the gas is a dilute Bose gas. 

Using the low-density (or small gas parameter) expansion allows the computation 
of the molecular chemical potential /ib, the speed of sound wave propagation u, the 
condensate depletion, and the ground state energy Eq. In the lowest order these are 

47rnbab ,„ ^ . 

^b = , (3.1) 

mb 



/47rnbab , . 



m^ 



no,h = f^b 



1 8^/nba^ 



3 
3 V vr 

2 



(3.3) 



Eo/V = '^^^. (3.4) 

rrih 

rrib = 2r7T, is the mass of a molecule and V is the volume of the system. The molecule- 
molecule scattering length ay^ appearing in these equations is proportional to the atom- 
atom scattering length. The coefficient of proportionality was first calculated by Petrov 
et al [59] who found 

ab « 0.60a. (3.5) 

This result was later confirmed in Refs. [14, 50] with the calculations of Ref. [50] 
presented in section 3.6 below. It should be noted that the gap equation predics Ob = 2a. 
Hence it is seen how the gap equation does not take interactions between the molecules 
properly into account. 
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At the crossover between BEC and BCS physics, the scattering length a diverges. 
This means that the gas parameter is no longer small, instead the scattering length 
will be much larger than the inter-molecular distance. The point at which a diverges is 
known as the unitary point. The physical behavior of the superfluid at zero temperature 
can only depend on one parameter, the density. The physical properties of the gas in 
this regime may be obtained by methods such as a large A^ expansion [73] , Monte Carlo 
techniques (see e.g. the review [32] for references), or other methods. The perturbative 
approximation employed in this thesis breaks down close to the resonance and here only 
systems with n}'^a <^ 1 will be considered. 

In this chapter, wide resonances are investigated. As discussed in chapter 2 this 
requires a ^ tq. It is important to note that while a diverges across a Feshbach 
resonance the effective range stays roughly constant. Thus the validity of the low- 
density expansion is restricted to scattering lengths such that |ro| <^ a <^ n^^'^. This 
limitation of the model's validity is not too restrictive; typical experimental densities 
of two-component Fermi gases are of order 10^^ cm^"^ corresponding to an interparticle 
distance of order lO^ao (ao is the Bohr radius) while r^ is typically of order SOoq. 

The fact that the investigated Feshbach resonances are wide implies that the 
dynamics of the closed-channel molecule is suppressed. In the limit of a vanishing 
effective range it is then simpler to describe the gas by a one-channel model rather than 
the two-channel model described in chapter 2. The equivalence of these two approaches 
will be demonstrated below. 

It should be reiterated that all calculations presented in this thesis proceed at 
zero temperature. A finite temperature extension of these results will not be discussed 
here. For further references see Ref. [2] and references therein. 

The systematic perturbative expansion in the gas parameter may be continued 
below the lowest order. The first two corrections to the results (3.1) and (3.4) depend 
only on the gas parameter after which further details of the interaction become impor- 
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tant. The discussion of these higher order contributions to the chemical potential and 
ground state energy will be studied in chapter 4. 

Another natural extension of the work presented here is to consider a two-component 
gas with a mass-imbalance. This is treated in chapter 5. 

This chapter is organized as follows. First, the one-channel model is introduced 
and the diagrammatic approach to the BEC regime of the weakly bound molecules 
is discussed. Next, in section 3.4, physical properties of the BEC are computed in the 
lowest order gas parameter expansion. Finally, sections 3.5 and 3.6 contain the few-body 
diagrammatic calculations relevant for the BEC regime. 

3.1 The one-channel model 

In the previous chapter, the two-channel model was introduced to describe a sys- 
tem where a bound state of two fermions appeared as the magnetic field was tuned across 
a resonance. The Feshbach resonances described in this thesis (and indeed relevant to 
most BCS-BEC experiments) are wide, and for wide resonances it is more convenient 
to treat the gas of weakly bound molecules using what is known as the one-channel 
model. Below, it will be shown that these two approaches are equivalent in the limit of 
infinitely wide resonances. 

The one-channel model may be expressed by a functional integral 

Z = jvij^Vi^^Vij^Vij^e'^^. (3.6) 

The action S'f = 5*0 -|- S'int consists of a free part and an interacting part. The free part 
is 

+ ^U<x, (3.7) 



5o= E/'^'^^^^'^(4 + 2LI^ 



while the interacting part is 
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Here A is chosen positive because the interaction potential is attractive. The interactions 
described by S'int take place at one point in space which by itself is unphysical. The 
results should be independent of the precise short distance physics and thus S'int needs 
to be supplemented by a cut-off A. Note that the one-channel model described by Eqs. 
(3.6), (3.7), and (3.8) does not explicitly contain the closed channel molecule associated 
with a Feshbach resonance. 

3.1.1 Diagrammatic approach 

Consider now the one-channel model of a two-component Fermi gas with m^ = 
m^ = m and 11-^=11^= fi. The propagator of fermions in the model described by Eqs. 
(3.6), (3.7), and (3.8) is 

Goip, ^) = 2/0^^ ^-n ' (3-9) 

to — p'^/zm + fi + zU 

where m is the mass of either one of the fermionic atoms. 

The scattering length of fermions in the vacuum may be computed by summing 

the diagrams shown in Fig. 3.1. The result of this summation is the T-matrix and the 

relationship between the scattering length and the T-matrix is^ 

a = l^r(O) (3.10) 

with rrir = m/2 the reduced mass of two atoms. 

In the calculation of the fermionic scattering length it is important to keep ^u = 
since the scattering proceeds in the vacuum. The diagrams in Fig. 3.1 form a geometric 
series since each loop is independent of the other loops. The T-matrix evaluated at zero 
incoming and outgoing fermion four-momenta is 

T(0) = A + An(0)A + An(0)An(0)A + . . . = Y^^^- (3.11) 

^ The T-matrix is a quantity defined for any kind of scattering. In tliis tliesis, scattering between 
two fermions, a fermion and a bosonic molecule, and between two molecules will be considered. These 
result in T- matrices T , T^ , and T^^. In general the superscripts will be suppressed to ease notation 
except where any confusion might occur. 
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+... 



Figure 3.1: The diagrams which sum to give the atom-atom scattering matrix T. The 
straight hne is a propagator of fermions while the square at each vertex is the coupling 
A. 



The bubble diagram 11 consists of two fermion propagators integrated over the free 
four-momentum in the loop. That is, 



n(o) = ^y"^Go(p)Go(-p) 



7TT-A 



(3.12) 



Thus, in the one-channel model the atom-atom scattering length is 



Att 2A\~^ 
m\ vr ) 



(3.13) 



In the one-channel model, the dynamics of the closed-channel molecule is sup- 
pressed. Instead the propagator of the closed- channel model is constant, equal to —A. 
However, the molecule needs to be "dressed" by the self energy correction due to re- 
peated interactions between the two fermionic components. Fig. 3.2 shows the propa- 
gator dressed by the self energy corrections, which consist of fermionic bubble diagrams. 
The molecular propagator is then 

-A 



Do{p) 



1 + An(p) ' 



(3.14) 




+... 



Figure 3.2: The relation between the dressed molecular propagator (thick wavy line) 
and the bare molecular propagator (thin wavy line). The bare propagator equals —A. 
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where it has been used that the diagrams again form a geometric series. The calculation 
now proceeds in the medium at a finite chemical potential /U < 0. Then the fermionic 
loop takes the value 

'4, 



n(p,^) = ^|^G{p + q)G{-q) 



= --^ + ^— V-w +p2/4"i - 2^ - iO. (3.15) 

27r^ 47r 

The infinitesimal is derived from the fermionic propagators and is such that also the 
molecular propagator is retarded. Note the importance of having /U < 0. Had this 
not been the case, the fermionic propagator would not have been retarded and hole 
propagation would have occured. 

Using the relation (3.13) between the scattering length and the cut-off, the molec- 
ular propagator becomes 

47r 1 , , 

Do(p,a;) = , (3.16) 

m a-^ - ^/my'-uj + p^/Am - 2;U - iO 

3.2 BCS-BEC gap equation 

In order to derive the BCS-BEC gap equation (the equation for the extremum 
of the action) from the one-channel model, it is convenient to introduce a Hubbard- 
Stratonovich field A. The one-channel functional integral (3.6) may then be written 
as 

Z= I Vij^Vij^VitiVi^iVAVAe''^', (3.17) 

with the action S[^t replaced by 

'1 



Sns = - / d-^x dt 



^■AA + (AV^tV^i + AV'iV'T) 



(3.18) 



The advantage of introducing the Hubbard-Stratonovich field is that the action 
S'f becomes quadratic in fermion fields, and these may then be integrated out. The 
resulting functional integral is 

Z= fvAVAe''^^, (3.19) 
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with 



UJ 



Sa = —i tr log 



2m + /^ 

-A 



-A 



^^ + 1^-^' 



A 



d-^xdtAA. 



(3.20) 



Here lo~^ = uj + iOsign uj. 

The gap equation is then derived by assuming that the extremum of the action 
Sa occurs at a constant value Aq of A. This results in 

1 



1 _i r d^p 

\~2] (2^ 



(3.21) 



2m ^ 



+ A2 



Using Eq. (3.13) the divergence in the momentum integral of Eq. (3.21) at large 
momenta may be regularized, and the gap equation becomes 



m 
Ana 



d^p 

(2^ 



1 



I 2m 



Ml +A2 



2m 
p2 



(3.22) 



3.3 Equivalence of the one- and two-channel models 

The two-channel model described in section 2.2 contains the dynamics of the 
molecular field. Compared with the one-channel model, this model has an additional 
parameter, the detuning. This parameter allows for the control of the effective range tq 
independently of the scattering length a. The precise relationships are listed here again 
for convenience 



mg 



ro 



4- i^o - ^ 



mg^ 



(3.23) 

(3.24) 



To ensure that a Feshbach resonance is wide, the effective range should be vanishingly 
small. This means that the interaction strength must be taken to oo, while at the same 
time the detuning should be adjusted in such a way that the scattering length remains 
finite. 
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This may be achieved as fohows. Redefine the bosonic fields of the two-channel 
model 

A = gb, A = gb. (3.25) 

Using this definition, the interconversion part of the action, Eq. (2.14), becomes 

5bf ^ - / d^xdt [AV^tV'i + AV'i^t] • (3-26) 

For the free bosonic action use 

eo = -g [-, TT^ = — T' (3-27) 



^4™ 27rV A' 

which follows from Eqs. (3.13) and (3.23). Then in the limit g ^ oo the free bosonic 
action (2.13) becomes 

S'ob = - d^xdt AA, (3.28) 

and the sum S'ob + 'S'bf exactly matches the Hubbard-Stratonovich action Sus of the one- 
channel model, see Eq. (3.18). Thus in the limit of wide resonances, the one-channel 
model is equivalent to the two-channel model. 

In order to ensure that the Feshbach resonance is wide, the interaction strength 
is taken to the limit g ^ oo. Thus it is again seen that it is not possible to treat g as a 
small coupling and construct a perturbative expansion in powers of g. 

3.4 The BEC regime of weakly bound molecules 

Having set up the perturbative framework in terms of the gas parameter, the 
calculation of properties of the BEC of weakly bound bosonic molecules proceeds in a 
fairly standard manner [1, 29]. 

It is convenient to trade the expansion in powers of na'^ with a diagrammatic 
expansion in the number of propagators beginning or ending in the condensate. Each 
such line is assigned the value Aq, the expectation value of the field A. Ag is not the 
density of particles in the condensate (the value assigned to condensate lines in the 
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standard dilute Bose gas) but is proportional to it, with the relation given below. Since 
the density of bosons in the condensate is lower than the total density, the perturbative 
expansion in powers of the gas parameter may be viewed as an expansion in the number 
of condensate lines in the diagrammatic approach. 

The presence of the condensate leads to the appearance of both normal (one 
particle enters, one particle exits) and anomalous (two particles enter and no particles 
exit or vice versa) propagators. These replace the "bare" propagators Go and Dq and 
will be denoted G^, Gg, for the fermions and D^, Dg. for the bosons. The fermionic 
propagators are needed in order to properly compute the particle number while the 
molecular propagators are needed to determine parameters of the gas. In the systematic 
perturbative treatment, the low-density approximations to the normal and anomalous 
propagators must be constructed. 

To this end it is convenient to construct the self energy corrections to the propaga- 
tors. The normal and anomalous fermionic self energies are denoted Sn and Sa while Sn 
and S20 are the normal and anomalous self energies for the molecules. The perturbative 
expansion in the number of condensate lines then consists in letting S = S^-*^) -|- S'^' -|- . . . 
where the superscript denotes the number of condensate lines. However, as discovered 
by Lee, Huang, and Yang [45, 46] non-trivial orders may be obtained through the sum- 
mation of certain classes of infinite numbers of diagrams. This does not affect the lowest 
order results presented in this chapter and is instead the subject of chapter 4. 

No normal self energy for the fermions has only one condensate line and thus 

S(,^) = 0. (3.29) 

The anomalous self energy in this order is given by the diagram shown in Fig. 3.3. This 
has the value 

41) = Ao. (3.30) 

In the next order, the contributions to the normal self energy are as shown in Fig. 3.4. 
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Figure 3.3: Lowest order contribution in the low density approximation to the anomalous 
fermionic self energy, Sa- Straight lines symbolize the fermionic atoms while the dashed 
line is a condensate line. The external fermion lines do not contribute to the self energy, 
they are present only to guide the eye. 



The square block shown in this diagram equals the sum of all diagrams contributing to 
scattering of an atom and a molecule. There are an infinite number of such diagrams. 
The computation of such a sum will be discussed in detail in section 3.5 below. Denote 
by t^^ (the atom-molecule scattering matrix) the result of the summation.^ Then 

4^) = Agt^f. (3.31) 

No diagrams with two condensate lines contribute to the anomalous self energy and 
thus 

Sl^) = 0. (3.32) 

For the molecular self energies, no diagrams with only one condensate line con- 
tribute and therefore 

S;i) = 0, ^i]l = 0. (3.33) 

In the next order, the self energies are shown in Fig. 3.5. The square block is the 
t-matrix of molecule-molecule scattering calculated in section 3.6. Again this consists 
of an infinite number of diagrams which may be summed. The self energies are 

T,f^ = 2A^t^^ (3.34) 

e(? = Ay\ (3.35) 



^ The simplest such diagram consists of Ei repeated twice and should not be included (this is the 
diagram depicted in Fig. 3.6, excluding external lines). This technical point will not be important for 
the following discussion. 
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Figure 3.4: Lowest non-vanishing contribution to the normal fermionic self energy, Sq. 
The square block is the sum of all diagrams contributing to atom-molecule scattering. 



The factor 2 is due to the symmetry of the diagram. 

Using the bosonic self energies Sn and S20, the normal and anomalous propa- 
gators of the bosons may be computed. This was first done by Beliaev [8, 9] with the 
result 



Ai(p,u;) 
The denominator has the form 



-Cq \-P, -^) - Sll(-p, -UJ) 



Dip,u;) 



D{p,u;) 



(3.36) 
(3.37) 



D{p,uj) = -['E2o{p,uj)f + {Dq^{p,uj) -Sii(p,u;)) (Dq ^(-p, -w) - Sii(-p, -w)) . 

(3.38) 
In the BEC regime, it was shown by Hugenholtz and Pines [38] that the propaga- 
tor of bosons must have a pole as the four-momentum tends to zero, corresponding to 
the presence of a gap-less sound mode in the condensate. This is the crucial observation 
on which the results of this section depend, and in general the Hugenholtz-Pines relation 
takes the form 

L>(0,0)=0. (3.39) 

In the lowest order in the expansion in powers of Aq, the propagator of bosons reduces 
to the form Dq obtained in Eq. (3.16). The chemical potential in the lowest order is 
then 

^^ = --^. (3.40) 



2raa? 




29 



Figure 3.5: Normal (top) and anomalous (bottom) self energy contributions to the 
bosonic propagator in the lowest non- vanishing order. Wavy lines are molecular propa- 
gators. The square blocks are the sum of all diagrams participating in molecule-molecule 
scattering. As with the fermion self energies, the external molecular propagators are 
not a part of the self energies. On the other hand, the condensate lines do contribute 
to the self energies. 



Thus, in the lowest order of the BEC regime the chemical potential of fermions equals 
half the binding energy of the molecular state as expected. 

It should be noted that the same result for the chemical potential would have 
been obtained if in the gap equation (3.22) the expectation value Aq had been ignored 
compared with fi. 

The bosonic self energies should be related to the scattering length of the molecules 
in vacuum, a^- The self energies appearing in the Hugenholtz-Pines relation are on the 
other hand proportional to the molecule-molecule scattering t-matrix in the presence of 
the medium. In particular the t-matrix is to be evaluated at a finite chemical poten- 
tial /i. However, as long as corrections to the chemical potential are small, t in the 
medium will approximately equal t evaluated in vacuum and will thus (as the four- 
momentum tends to zero) be proportional to the vacuum scattering length Ob. The 
precise relationship is 



.bb 



27r 
mZ 



ah- 



(3.41) 
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Figure 3.6: The diagram which leads to the lowest order particle number (3.47). As 
opposed to the self energy diagrams, here the external lines must be included. 



Here 

Z = ^ (3.42) 

is the residue of the molecular propagator at its pole and is needed for proper normal- 
ization of external propagators [10]. 

The chemical potential of the bosonic molecules is defined as 

/ib = 2/i+ 2' (3-43) 

which vanishes in the lowest order of the low-density approximation. Using this def- 
inition of the bosonic chemical potential, the Hugenholtz-Pines relation (3.39) results 
in 

/xb = Z [Sn(0) - S2o(0)] - ^^^ [Sn(0) - S2o(0)]' . (3.44) 

Eq. (3.39) will have one additional solution, however this solution is unphysical [1, 29]. 
In the lowest order in density 



^b ~ Z[Sn(0)-S2o(0)] 
aahmAQ 



(3.45) 



and the relation becomes reminiscent of the Hugenholtz-Pines equation in the standard 
dilute Bose gas [38]. 

The particle number may now be computed from the fermionic self energies. The 
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particle number equation is 



= -^ 1^^ I ^[Gn,^iP,-^) + Gn,i{P,^)], (3.46) 

where C is the contour which includes the real axis and is closed in the upper half 
plane ensuring that the field operators occur in the proper order. The diagram which 
contributes to the particle number in the lowest order is shown in Fig. 3.6. Thus the 
particle number is 

= — -. — -■ (3.47 

4tt 

The sign of the expression comes from fermionic anti-commutivity and the factor 2 from 

the two species of fermions. In the evaluation of the integral, the chemical potential has 

been set to its zero'th order value, Eq. (3.40). With the help of Eq. (3.47) the bosonic 

chemical potential Eq. (3.45) becomes 

vrabn 47rabnb ,„ ,„. 

/Ub = = • (3.48) 

m nib 

The ground state energy Eq (measured from the binding energy of the molecules, 

i.e. from —nij/ma?) is calculated using the standard relation [38, 29] 

Eq/V - ^nb/ib = J^l^J (^ (w +//4m) D^{p,u), (3.49) 

where the contour C is again taken in the upper half plane. In the lowest order, the 
normal molecular propagator reduces to Dq which is a retarded propagator. Thus with 
the help of Eq. (3.48) the ground state energy becomes 

^o/^=^nbMb = ^^^^. (3.50) 

2 m\) 

The effect of including the higher order molecular propagator (3.36) in the equation for 

the ground state energy will be discussed in chapter 4. 
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Even at zero temperature there is still a finite fraction of the bosonic molecules 
not in the condensate. This phenomenon, caused by the presence of interactions, is 
known as the condensate depletion. The total density of bosons is 

duj I' d?p 
Z Jc2tt J J2n) 



nh = no^b + — I 7:^ I 7w:^Dri{p,^), (3.51) 



with riQ^b the density of bosons in the condensate. In the lowest order, the density 
'^b = '^cb- Using Eq. (3.36) for D^i(p,uj) it is seen that the frequency and momentum 
dependence of the molecular self energies T,ii(p,uj) and T,2o{p,uj) is needed. However, 
considering the integral in Eq. (3.51), it is observed that the integration is dominated 
by frequencies of order cj ~ ^Ub ~ ^0 ^^'^ momenta of order p^ /Am ~ Aq. For such low 
frequencies and momenta, the underlying structure of the molecular propagator (3.16) 
is not important and it reduces to 

^o(p, ^) « 277^T ^l,- (3-52) 

uj — p'^ /Am + ;Ub + *0 

Also, ZSii(O) = 2/ib and ZS2o(0) = /Ub coincide with Sii(O) and S2o(0) in the usual 
dilute Bose gas. Hence, for small frequencies and momenta, the normal molecular 
propagator divided by Z coincides with the dilute Bose gas normal propagator at this 
order in the expansion. The result of evaluating the integral (3.51) then takes the usual 
form [1, 29] 

no,b = nb 1--W^^ . (3.53) 

o y vr 

Finally, the speed of sound wave propagation is evaluated by use of the Hugenholtz- 
Pines relation, Eq. (3.39), at finite p,LO. For small uj < ^\^ and p^/Am < ^Ub the self 
energy terms approximately take the limiting values S^-^ (0) and S20 (0) in the low 
density expansion. Imposing the condition 

D{p,Lo)={) (3.54) 

for the denominator, Eq. (3.38), of the molecular propagators, it is found that 



2 2 2 /4vrnbab , . 

Lo = u p , u=i K — , (3.55) 

V "^b 
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where u is the speed of sound. 

The results of this section coincide with the results for the dilute interacting Bose 
gas [1, 29] with scattering length Cb- In order to go beyond these results it is necessary 
to modify the self energies in a systematic way. This will be further investigated in 
chapter 4. 

3.5 Scattering of a molecule and an atom 

In this section, the scattering problem of a fermionic atom with a bound state 
of fermions is treated. The solution of this scattering process has a long history. It 
was first solved in the approximation of vanishing effective range by Skorniakov and 
Ter-Martirosian in 1956 [68] who studied the related problem of a neutron scattering 
with a deuteron (a neutron-proton bound state). In the context of cold two-component 
Fermi gases the scattering problem was solved in Ref. [56] for any mass ratio of the 
fermionic atoms below 13.6. The three-body calculation presented below is a useful 
exercise before proceeding to more complicated few-body calculations. It is important 
to note that the calculation below relies on the approximation of a vanishing effective 
range, i.e. a ^ jro] and the results are only valid in this regime. 

It should be mentioned that there is a fundamental difference between the fermionic 
and the bosonic problem. The system of three identical bosons with a large scattering 
length has an infinite number of bound states, the effect known as the Efimov effect 
[25]. In this case, to correctly describe the three-body physics, an additional three-body 
parameter coming from short distance physics is required [20, 7, 12]. This is not the 
case in the fermionic problem (for mass ratios below 13.6); here the three-body problem 
is completely describable in terms of the two-body scattering length and the fermionic 
masses [56]. 

In the one channel model, the only length scale left is the scattering length between 
the distinguishable atoms. Since no additional parameters are required in order to 
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Figure 3.7: The Born approximation to atom-molecule scattering. 



describe the three-body problem, the atom-molecule scattering length, Obf, must be 
proportional to the atom-atom scattering length, a. 

Apart from the three-body scattering length, it is also interesting to note that 
the achievable densities for two-component Fermi gases are limited by three-body re- 
combination, a process governed by three-body physics. Here two atoms form a deeply 
bound state with binding energy of order l/mR"^ and a third atom escapes the trap 
as it carries away the released binding energy as kinetic energy. The process requires 
at least two identical fermionic atoms to approach each other closely (at the order of 
the range of short distance physics. Re) and is suppressed by the Pauli principle. The 
loss rate depends on the precise three-body physics and may be extracted from the 
three-body calculations below. A precise discussion will be postponed for chapter 5 on 
mass-imbalanced two-component Fermi gases since the loss-rate depends intimately on 
the mass ratio. 

This is the first of a number of few body problems treated using integral equations 
and therefore this will be done very carefully. The simplest process contributing to the 
atom-molecule scattering is that the molecule breaks into its two constituent atoms, one 
of which forms a bound state with the remaining atom. The process is illustrated in 
Figure 3.7; had the scattering been weak, this would have been the only contribution 
to the scattering and thus this diagram will be called the Born contribution. However, 
as will be discussed below, all of the diagrams contributing to the scattering are of the 
ladder diagram type shown in Fig. 3.8 and these are all of the same order. 

As pointed out above, the vacuum scattering diagrams are exactly those needed 
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Figure 3.8: The ladder type diagrams contributing to atom-molecule scattering. 



for the normal fermionic self energy at next to leading order in the low density expansion 
and the normal and anomalous molecular self energy at lowest order. That the type of 
diagrams shown in Fig. 3.8 is the only type contributing to the scattering stems precisely 
from the fact that this is vacuum scattering, i.e. all propagators are retarded and there is 
no hole propagation. This in turn means diagrammatically that all propagators specify 
a time direction and the possible diagrams are greatly simplified. In fact, because of 
the simplicity of the diagrams it is possible to sum all the possible diagrams by using 
an integral equation of Lippmann-Sch winger type. 

Since the scattering proceeds in the vacuum, the chemical potential n will be 
taken to vanish in this section. 

The scattering t-matrix consists of all diagrams with external incoming and out- 
going lines being the scattering particles. The external lines are not included in the 
t-matrix, however they do affect the overall normalization. The first diagram in Fig. 
3.8 thus contributes a single atomic propagator, proportional to ma^. In general, a 
diagram with n loops will then contain 2n + 1 fermionic lines of order (ma^)^"^"'^, n 
bosonic lines of order (a/m)'^, and n integrations over four-momenta ~ rTi~"a~^". It 
follows that each diagram will be of order ma'^. 

Given the structure of the contributing diagrams it is possible to write an integral 
equation for the sum of these. The integral equation is shown in Fig. 3.9. In Eq. (3.13) 
the two-fermion bubble diagrams were summed to yield the atom-atom scattering length 
by using a geometric series. The reason why the three-body problem is an integral 
equation rather than a geomatric series is that there will be an integral over the free 
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Figure 3.9: The integral equation for the t-matrix (shown as the square block) of atom- 
molecule scattering. 



four-momentum in the loop of the iterated term. 

Let the four-momenta be as shown in Fig. 3.9. That is, the incoming molecule 
and atom have four-momenta {p,po + E) and {—p, —po)i respectively, while the outgoing 
molecule and atom have four-momenta {q, qq + E) and {—q, —qo), respectively. E is the 
total energy going into the scattering and since the initial and final states consist of a 
bound pair and an atom, it is necessary to restrict the total energy to energies E < E]^. 
Otherwise, the final scattering state could be three free particles. The scattering t- 
matrix with these kinematics is denoted t{p,pQ;q,qQ) (the dependence on total energy 
is suppressed in this notation). Then the integral equation becomes 



t{p,po; q, qo) = -G{p + q + E)-i 



(27r)4 



G{Q + q + E)G{-Q)D{Q + E)t{p,po; Q, Qo). 

(3.56) 



Here, in a slightly abusive notation, the total energy E is also used as the total four- 
momentum {0,E). To simplify notation further, the four-momentum {p,po) is abbrevi- 
ated as p. The propagator of molecules is here given by the expression (3.16) with the 
chemical potential taken to vanish. 

The t-matrix depends on five variables; the two frequencies, the amplitudes of the 
two momenta, and the angle between incoming and outgoing momenta. However, as 
the dependence of t on p and po is the same on both sides of the integral equation, the 
integral equation for t is only a three dimensional problem. Below, it will be discussed 
how the frequencies qo and Qo rnay be integrated out from Eq. (3.56) and in chapter 5 
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it will be demonstrated how the integral equation may be projected onto each spherical 
wave component, making the solution of each of these a one- dimensional problem. 

^(P)PoJ 9) Qo) is analytic in the upper half planes of both pQ and go- This may be 
seen from the ladder diagrams contributing to the t-matrix or simply from Eq. (3.56). 
This in turn means that the integration over Qo in Eq. (3.56) may be carried out by 
closing the contour in the upper half plane, with the only contribution coming from the 
simple pole of G{—Q), setting Qo = —Q'^/2m. The solution of the integral equation is 
then most easily obtained by letting qq -^ —q^ /2m which means that the momentum 
dependence of the t-matrix is the same on both sides of the integral equation. Should 
the value of the t-matrix be desired at a different value of gO) the result of solving the 
integral equation at qq = —q^ /2m may simply be inserted on the right hand side of Eq. 
(3.56) and the integration carried out. In this sense the t-matrix with po = —p^ /2m, 
and qq = —(^ /2m is quite natural. It is the result of attaching the external atomic 
propagators to the t-matrix and integrating over their frequencies, closing the contour 
in the upper half plane (where t is analytic) . This will be called the "on-shell" t-matrix 
and satisfies the integral equation 

t{'p,-p^/2m;q,-q^/2m) = -G{p + q + E) 

- j ^%G{Q + q + E)D{Q + E)t{p, -pV2m; Q, -Q'' /2m). (3.57) 

The iterative procedure described above for obtaining the solution to the integral equa- 
tion at the desired values of pO)9o from the on-shell t-matrix is illustrated in Fig. 3.10. 
To calculate the s-wave scattering length, the total energy should be set to the 
two-body binding energy. The scattering length is then found by solving the integral 
equation, averaging over angles, and taking the limit of vanishing incoming and outgoing 
four-momenta. Equivalently, the incoming momentum can be taken to vanish prior to 
solving the integral equation after which there can be no dependence on the angle 
between incoming and outgoing momentum. Using this second approach yields the 
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Figure 3.10: The procedure for relating the on-shell t-matrix (the square block on the 
right hand side) to the t-matrix evaluated at a desired frequency (the left hand side). 



integral equation 



tip) =^+^r d, ^^('^^ log '+pi+^i+<ip. 

'■'^' l+p2 TrpJo ^l-y/l + 3Q2/4 ^1+p^+q^-qp 



(3.58) 



In this equation, momenta are measured in units of inverse scattering length and the 
"on-shell" s-wave scattering matrix is t(p) = t(0, 0;p, —p^ /2m). 

In general, the integral equations have a form which in matrix notation is 



T'l — Vi -\- l\ijTj. 



ir3- 



(3.59) 



The matrix equation has the formal solution 



U = {I-K)T^\,. 



(3.60) 



The matrix / is the identity matrix. The numerical solution of the t-matrix is then 
a matter of choosing a discretization and solving the matrix equation. Typically, for 
purposes of convergence, it is advantageous to perform a change of variable such that 
the limits of integration become finite [61]. A convenient choice is 

2 



p 



z + l 



1, ^G]-l,l]. 



(3.61) 



All integrals are then performed employing Gauss-Legendre quadrature [61] . Using this 
method the t-matrix is found to take the shape shown in Fig. 3.11. 
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Figure 3.11: The t-matrix t(p) in units of ma'^ calculated by solving Eq. (3.58). For 
p <^ a~^ the t-matrix is approximately constant while for p ^ a^^ the t-matrix satisfies 
a power-law behavior, decaying as p~3.i7 ^jj_ 



The atom-molecule scattering length Obf is proportional to the t-matrix of Eq. 
(3.58) evaluated at vanishing momentum and energy E = E]^. To obtain the scattering 
amplitude from the t-matrix it is necessary to renormalize the external propagators. The 
fermionic propagators are already free propagators, however the molecular propagators 
need to be renormalized by the square root of the residue of the molecular propagator 
at the energy of the bound state [10], i.e. 



T(0) = Zt(0), 



(3.62) 



with Z given by Eq. (3.42). The scattering amplitude is related to the scattering length 

by 

(3.63) 



27r 3-K 

T[0) = Obf = — abf, 

m3i- m 



where m^^ is the three-body reduced mass. Using the method described above it is 
found that 

Obf ~ 1.18a, (3.64) 

in perfect agreement with the literature [68, 56]. 
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Bound states in the three-body problem show up as poles of the scattering am- 
plitude. This corresponds to solutions of the homogenous integral equation 

ti = Kijtj, (3.65) 

or in other words to the integration kernel K having a unit eigenvalue. Bound three- 
body states will have E < E\,, and writing E = (1 + e3)E\, the homogenous equation 
becomes 

tip) =^r d, ^ ^^(^^ log ^+^^+^^^;+^^ (3.66) 

vrpio 1 - Vl + £3 + 3g2/4 ^1 + es+p^+q^-qp 

It is found that the three-body problem does not have any s-wave bound states. Note 
that the validity of this statement is limited by the zero effective range approximation. 
There may still be deeply bound states with binding energies of order 1/mrQ. 

3.6 Molecule-molecule scattering 

As discussed above, properties of the gas of weakly bound molecules depend 
intimately on the molecule- molecule scattering length, Ob. Mean field theory gives the 
result Ob = 2a. This result implicitly assumes that the coupling is weak and that only the 
diagram shown in Fig. 3.12 contributes to the scattering. An approximate diagrammatic 
method was developed by Fieri and Strinati in Ref. [60] who found Ob ~ 0.75a. Later 
the problem was solved exactly by Fetrov, Salomon, and Shlyapnikov [59] by solving 
the quantum mechanical four-body problem using the Bethe-Feierl's method. In this 
work it was found that Cb ~ 0.6a. This last result is confirmed below by the use of 
diagrammatic methods. While the work presented in this section and published in Ref. 
[50] was in progress, the problem was solved by Brodsky et al in Ref. [14] by essentially 
the same diagrammatic method. 

The scattering t-matrix consists of all possible diagrams with two incoming and 
two outgoing molecules, excluding the external propagators. In the atom-molecule 
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Figure 3.12: The simplest contribution to molecule-molecule scattering. 



scattering problem treated above, it was discussed how diagrams contributing to the 
molecule-molecule scattering do not contain any internal condensate lines nor any hole 
propagation. This is again the case, since the process considered is scattering in the 
vacuum. This observation greatly simplifies the diagrams possible in the scattering pro- 
cess. Since all propagators can be assigned a definite time direction, at any time in the 
scattering of two molecules there will be either two molecules, one molecule and two 
atoms, or four atoms. 

With this observation in mind, consider the different diagrams which contribute 
to the scattering t-matrix. The simplest possible diagram is the Born diagram shown in 
Fig. 3.12 where the two molecules exchange their constituent fermionic atoms. Other 
possible diagrams contributing to the scattering are shown in Fig. 3.13. Fig. 3.13a shows 
the propagation of an intermediate molecule. Fig. 3.13b a slightly more complicated 
diagram, and 3.13c a diagram where a pair of molecules propagate at an intermediate 
step. The crucial observation is that all of these diagrams are of the same order, namely 
they are all proportinal to rrfia^. Indeed this is true for all possible vacuum scatter- 
ing diagrams. Thus the problem of molecule-molecule scattering is not amenable to a 
perturbative treatment. Instead of considering only the Born diagram. Fig. 3.12, an 
infinite number of diagrams need to be taken into account. 

Fortunately it is possible to perform the summation of all diagrams by using 
an integral equation. To this end, it is useful to define the concept of a two-boson 
irreducible diagram as a diagram which cannot be cut in two by cutting two molecular 



(a) 



(b) 
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(c) 
Figure 3.13: Diagrams which contribute to molecule-molecule scattering. 



lines only. For instance, the diagrams in Fig. 3.13a and b are two-boson irreducible 
while the diagram in Fig. 3.13c is not. The scattering t-matrix may then be viewed 
as the repeated propagation and scattering of two molecules, each pair of molecular 
propagators separated by the sum of all two boson irreducible diagrams. In this sense, 
the sum of two-boson irreducible diagrams is similar to the bubble diagram, 11 (see Fig. 
3.1) occurring in the repeated atom-atom scattering. However, there will be an integral 
over a free four-momentum and thus the sum does not form a geometric series. Defining 
r as the sum of all two-boson irreducible diagrams, t satisfies an integral equation as 
shown in Fig. 3.14. 

Each diagram contributing to F will have a number, n, of molecular propagators, 
4 + 2n fermionic propagators, and n + 1 integrations. The total order of such a diagram is 



[^n-^-T'^'i^r'--'"' 



Tma" ) ( ;:::br ) = w, a . Any diagram contributing to t will then have n+1 



insertions of F, 2n molecular propagators and n integrations and will be proportional 
to m^a?. 
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Figure 3.14: Integral equation for the i-matrix (shown as the square block) of molecule- 
molecule scattering. Here, for simplicity E has been used as the four-momentum (0, E). 



The integral equation for the t-matrix with kinematics as shown in Fig. 3.14 is 

*(P, Po) = r(0, 0; p, Po) + ^ <fd<l dqoD{q, qo+E/2)D{q, -qo+E/2)T{q, qo; p, po)t{q, qo). 

(3.67) 
The incoming momentum has been taken to zero while the outgoing four-momenta of the 
molecules are (±p, E/2zizpo). This choice of kinematics is convenient since ultimately the 
quantity of interest is the s-wave scattering length, proportional to t(0,0). The scalar 
t{p,pQ) cannot depend on the direction of p, which is why it is possible to average the 
sum of two-boson irreducible diagrams over the angle between incoming and outgoing 
momentum. This average is denoted T{p,po;q,qQ). To determine the scattering length, 
let E = 2Eb, while to search for bound states look for solutions of the homogenous 
integral equation with E < — |2£'b|. 

Because of the symmetry of the problem, it is convenient to slightly redefine what 
is meant by G and D. Below, in this section only, the convention will be used that all 
atomic propagators have their frequency shifted by a quarter of the energy going into 
the scattering, while the molecules have their frequency shifted by half the energy. That 
is 



G{p,Po) 
D{p,Po) 



Po - _p2/2m + E/4 + iO ' 
47r 1 



(3.68) 
(3.69) 



in a-i - ^Jm^-po - E/2 + p'^/Am - iO ' 
For the purpose of calculating the molecular scattering length a\^, the total energy is 
2£'b and in this case each molecule will carry the molecular binding energy. 
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Figure 3.15: Integral equation for the vertex H. Here the notation E has been used to 
denote the four-momentum {{),E). 



The integral equation (3.67) is not sufficient to compute the t-matrix. The sum of 
all two-boson irreducible diagrams still needs to be computed and this sum may again 
be obtained through the solution of an integral equation. It is not possible to directly 
write down an integral equation satisfied by T. However, it is possible to relate T to 
the solution of an integral equation for the sum, H, of two-boson irreducible diagrams 
with two incoming molecules, one outgoing molecule and two outgoing fermions. The 
integral equation to be solved is depicted in Fig. 3.15. This integral equation is quite 
similar to the equation for atom-molecule scattering, see Fig. 3.9. 

Let the kinematics be as shown in Fig. 3.15. H with incoming four-momenta 
of {^q,E/2 lb go) for the molecules, {—p,E/2 — po) for the outgoing molecule, and 
(p/2 ibQ, E/A +pq/2 ibQo) for the outgoing fermions is denoted E{q;p/2 + Q;p/2 — Q). 
The integral equation for H may then be written as^ 



E{q;p/2 + Q;p/2-Q) 
d^Q' 



Att 



[G{-p/2 - q + Q)G{-p/2 + q-Q) + {q^ -q)] 



I Wf ^^^^^^ " ^'^^^"^^/^ + Q')D{-P -Q + Q')E{q;p/2 + Q;p/2 - Q') 

+G{p/2 + Q')G{-?,p/2 - Q')D{-p + Q- Q')E{q; p/2 + Q';p/2 - Q) }.(3.70) 



^ The calculation of the vertex H presented here differs slightly from the formalism in Rcf. [50] . The 
present formalism is better suited to an extension to treating the mass imbalanced problem investigated 
in chapter 5. 
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Figure 3.16: The relation between T and H. 



The integration over directions of g (with notation j dO,^) is the averaging over directions 
of incoming momentum, i.e. the projection onto the s-wave. Minus signs in front of 
both the Born terms and the iterated terms all follow from anti-commutivity of the 
fermionic atoms. 

Eq. (3.70) has to be solved at every \q\ and qq. It is an integral equation in 5 
variables, the length of p and Q, the angle between these vectors, and the corresponding 
frequencies. As will be shown below, it is possible to integrate out the two frequencies 
such that the resulting integral equation becomes three-dimensional. 

Finally, the two-boson irreducible vertex F is related to H by joining the two 
external fermionic propagators into a molecule as shown in Fig. 3.16. This relation may 
be written as 

r{q,qo;p,Po) = lj ^G{p/2 + Q)G{p/2 - Q)E{q;p/2 + Q;p/2 - Q). (3.71) 

The factor 2 ensures that all diagrams enter F with the correct weights. 

It is advantageous to split 'r.{q]p/2 + Q]p/2 — Q) into parts 'B'^(q;p/2 + Q;p/2 — Q) 
[E^{q;p/2 + Q;p/2 — Q)] analytic in the upper [lower] half planes of Qq, with 

E{q;p/2 + Q;p/2 - Q) = E+{q;p/2 + Q;p/2 - Q) + E-{q;p/2 + Q;p/2 - Q). (3.72) 

To this end, note that the iterated terms are already thus split. To split up the Born 
terms, use the relation 



G{p)G{q) 



G{p)+Giq) 

G(p)-i + G(g)-i' 



(3.73) 
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Then Eq. (3.71) becomes 



1 f d?Q 



(27r) 



G{p/2-Q)^-{q;p/2 + Q;p/2-Q)l 



)o=-po/2+(p/2+Q)2/2m-£/4 



+G{p/2 + Q) E+{q;p/2 + Q;p/2 - Q)\q^__ 



=po/2-(p/2-Q)2/2m+B/4 



. (3.74) 



In each of the two last terms of Eq. (3.70) it is possible to integrate out Q'q by 
closing around the single pole in the lower half plane. The two parts of H appearing in 
Eq. (3.74) satisfy the coupled integral equations 

E- (q-p/2 + Q, (p/2 + Qf/2m - E/A-p/2 -Q,po- (p/2 + Qf/2m + E/4 



d?Q' 



G{-p/2 + q-Q) 



+ {q^ -q) 



(2vr) 



G'\-p/2 -q + Q)+ G~\-p/2 + q-Q) 
3 G(-3p/2 + Q', -po - {p/2 - Q'f/2m + E/A)D{-p -Q + Q') 



xEi q-p/2 + Q, (p/2 + QY/2m - E/A-p/2 - Q\ {p/2 - Q'f /2m - E/A 



E+ q-, p/2 + Q,po- {p/2 - Qr/2m + ^/4;p/2 - Q, {p/2 - QY /2m - E/4 



4tt 
(27r)3 



G{-p/2 -q + Q) 



+ {q^ -q) 



G-\-p/2 -q + Q)+ G-\-p/2 + q-Q) 

G{-3p/2 - Q', -Po - {p/2 + Q'f /2m + E/A)D{-p + Q - Q') 



xE q-,p/2 + Q', {p/2 + Q'y/2m - E/4-,p/2 - Q, {p/2 - QY /2m -E/A). 



(3.75) 

In the first of these equations Qq = —po/2 + (p/2 + QY /2m — E/4 and Q'q = po/2 — 
{p/2 - Q'Y/2m + E/4. In the second equation Qo = po/2 - {p/2 - QY /2m + E/4 and 
Q'o = -po/2 + {p/2 + Q'Y/2m - E/4. 

On the right hand side of Eqs. (3.75), the vertex E only appears in the form 



S(9;Pi,Pi/2m - E/4-,p2,pl/2m - E/4) = xiQ;Pi,P2)- 



(3.76) 



This will be called the "on-shell" vertex. The reason for this terminology is clear; 
had the propagators G{pi) and G{p2) been attached to the vertex H(g;pi;p2) and the 
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frequencies been integrated over using the poles of the propagators, then these on-sheh 
values of the frequencies would have been obtained. In this sense, the on-shell vertex is 
equivalent to the atom-molecule scattering vertex t^^{p,—p'^ /2m;q,—q^ /2m) discussed 
above, see Eq. (3.57). 

Solving the integral equations is now possible with the substitution 

PQ -^ {p/2 + Qf/2m + {p/2 - Qf/2m - E/2, (3.77) 

at which value also the left hand sides of Eqs. (3.75) will be on-shell. Adding the 
expressions in Eqs. (3.75) at this value of pQ results in an integral equation in three 
variables for the total on-shell vertex which may subsequently be inserted on the right 
hand side of Eqs. (3.75) to find the value of the vertex H at any value of pq. 

After a change of variables to pi = p/2 + Q and p2 = p/2 — Q, the equation for 
the total on-shell vertex takes the form 



x{q;Pi,P2) = - I ^[G{q-pi)G{-q-p2) + {q^ -q)] 

f d^Q { 

- j j^-r^ { G{-pi -p2- Q)D{-Q - pi)x{q; Pi, Q) 



+ G{-pi-p2-Q)D{-Q-p2)x{q;Q,P2) ,,, ^ ,(3.78) 

where {pi)o = Pi/2m — E/A and (^2)0 = pI/^'^ ~ E/A. 

Inserting Eq. (3.78) into the right hand sides of Eqs. (3.75) and the resulting 
expression into the right hand side of Eq. (3.74), the vertex T finally becomes 

r(g,go;p,po) = r(°)(g,go;p,po) 

1 f d^pi cP'p2 



Here again {pi)o = p\/2m — E/A and (p2)o = Pi/2m — E/A while it is important to 
note that pQ is a free parameter. This method (3.79) of constructing the two-boson 
irreducible vertex corresponds diagrammatically to the situation depicted in Fig. 3.17. 
r*^^' is the result of calculating the Born diagram, shown in Fig. (3.12), 





Figure 3.17: An illustration of Eq. (3.79). 
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dJl^ 



d^Q 



G{Q + p/2 + q/2)G{Q - p/2 - q/2) 



At: J (27r)4 

xGi-Q - p/2 + q/2)G{-Q + p/2 - q/2) (3.80) 



d^Q 



A 



dr^n 



47r J (27r)3(^2_^2)(^2_C'2)^ 



A = E/2 - Q'^/m - p^/Am - q^/Am, 
B = Po-Q ■ q/m, 
C = qo-Q- p/in. 



(3.81) 



(3.82) 



The angular integrations may be performed independently by first integrating over 
directions of g and then over directions of Q. The result is 

^2 foo ^Q 



T'-°'>{q,qo;p,Po) 



m 



IGiT^pq Jq —E/2 + Q"^ /m + p^ /Am + q^ /Am 
X lo 



{-E/2 + Q'^/m + p^/Am + q^/Am + Qq/2mf - pi 



X loe 



{-E/2 + Q'^/m + p^/Am + q^/Am - Qq/2mY - p^ 
{-E/2 + Q'^/m + p'^/Am + q'^/Am + Qp/2m)'^ - ql 



(3.83) 



{-E/2 + Q'^/m + p^/Am + q^/Am - Qp/2mY - ^q 

The procedure for solving the set of integral equations is now as follows. The 
integral equation (3.78) for the on-shell x has to be solved at every {q^qo). However, 
the kernel of this integral equation does not depend on (g, go)) an observation which is 
of great help in the numerical solution, since this means that the identity matrix minus 
the integration kernel needs to be inverted only once [see Eqs. (3.59) and (3.60)]. With 
the on-shell x determined, this should be inserted in Eq. (3.79) to find the two-boson 
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Figure 3.18: The original and rotated contour of frequency integration used in Eq. 
(3.67). The crosses correspond to poles of the molecular propagators while the black 
lines correspond to the square root branch cuts. The gray boxes are the areas in which 
the two-boson irreducible vertex T may have non- analytic structure. 



irreducible vertex T. Finally, F may be inserted in Eq. (3.67) to find the scattering 
t-matrix. 

A bit of care needs to be taken in peforming the integration over frequency in 
Eq. (3.67). Due to the square-root structure of the molecular propagators there will be 
branch cuts starting at 



^0 



-E/2 + q^/Am > 0, and qo = E/2 - q^/Am < 



(3.84) 



and continuing towards ibcx), respectively. The contour of integration is such that for 
Re(go) > [Re(go) < 0] the integration is above [below] the branch cut, which follows 
from the infinitesimal in the molecular propagator, Eq. (3.16). Furthermore, the molec- 
ular propagators have poles at go = ±(g^/4?n, + E\y — E/2 — iO). It is thus advantageous 
to perform a Wick rotation as illustrated in Fig. 3.18. This removes the problem of 
integrating close to the branch cuts and poles. In order to perform such a Wick rotation, 
the vertex F(p,po; 9, ^Zo) needs to be analytic in quadrants / and ///. Eqs. (3.78, 3.79, 
3.80) show that this is indeed the case and that any non-analytic structure of F must 
be restricted to the areas marked as gray boxes in Fig. 3.18. 

When computing the scattering length, the total energy should be set to E = 2E\y. 
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Then the simple poles of the molecular propagators present in Eq. (3.67) occur at 
qo = ib(g^/4?72 — iO), i.e. close to both the original and the Wick rotated integration 
contour. These simple poles are integrable and may be treated by a change of variables 
to "polar" coordinates, 

q'^/4m = R^ sin 9, qq = R^ cos 6. (3.85) 

Here R G [0, oo[ and 9 G [0, tt]. 

To search for bound states in the four-fermion problem, let E < 2E]^ and look 
for solutions to the homogenous equation. In a matrix notation this corresponds to the 
existence of a solution to the equation ti = Kijtj or, in other words, the existence of a 
bound state implies that the kernel of the integral equation, K, has an eigenvalue equal 
to 1. No bound states are found by using this method. As the one-channel model is 
employed here, the effective range vq has been taken to zero. Thus the validity of this 
method is only for bound states with binding energies |-E| ^ — i and the result does 
not exclude deeply bound states. 

The scattering amplitude is calculated on-shell, with the incoming and outgoing 
molecules all having four-momentum (0,E}^). Thus E = 2E]^. As in the three-body 
problem above, to calculate the scattering amplitude each external molecular propagator 
needs to be renormalized by yZ, and the relation between the scattering amplitude and 
the scattering matrix is 

T(0) = Z2t(0,0). (3.86) 

The scattering length is then found through 

27r 
m = -a, (3.87) 

with Mr = m the molecule-molecule reduced mass. By solving the set of equations 
described above, it is found that 

Ob « 0.60a (3.88) 
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Figure 3.19: The scattering amplitude Z'^t{p,po) in unites of a/m evaluated at po = 0. 



in complete agreement with Refs. [14, 59]. Fig. 3.19 shows the computed t-matrix, 
from which t(0, 0) may be extracted. 



Chapter 4 
Quantum corrections 



In chapter 3 properties of the BEC regime of the dilute Bose gas consisting of 
weakly bound molecules were calculated. The calculations were performed in the low 
density limit and it was seen how this limit corresponds to a diagrammatic expansion 
in terms of lines beginning or ending in the condensate. The results obtained were 
seen to exactly match the results of the standard dilute Bose gas [1 , 29] with the effect 
of the fermionic atoms being that the boson-boson scattering length is related to the 
fermion-fermion scattering length by a^ ~ 0.60a. 

It is possible to extend the low density expansion in order to obtain the chemical 
potential and ground state energy in higher orders. This involves the determination of 
the fermionic and bosonic self energies in higher orders. Somewhat surprisingly, the next 
order self energies for the bosons do not involve diagrams such as in Fig. 4.1 which only 
contributes at even higher order. Rather, they arise due to non-trivial resummations 
of diagrams occuring because of infrared divergences of the molecular propagator. Lee, 
Huang, and Yang computed the first correction to the ground state energy [45, 46]. 
Later, results were obtained in the next order by Wu [76] and Sawada [65] and the 
results of these authors confirmed by Hugenholtz and Pines [38] . 

The results in the dilute Bose gas are the following: 



47rabnb 

Mb = 

"lb 



32 n^al A y/?,\ n . 



(4.1) 
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Figure 4.1: A high order contribution to the molecular self energy. 



Eq _ 2TTa\,nl 



V 



TTlb 



128 Rn,/ 4 

15 V ^ I 3 



^" 



■K 



ribflb log nbCb 



(4.2) 



Below, it will be shown that these results are also correct in the gas of weakly bound 
molecules, with the scattering length a\^ related to the underlying fermion scattering 
length a. 

The first two corrections to the ground state energy and chemical potential are 
universal in the sense that they only depend on the scattering length Ob, the density, and 
the mass of bosons. Higher order corrections will depend on the shape of the interatomic 
potential, in particular on the effective range tq. These will not be considered here. 

4.1 The molecular propagator at low densities and momenta 



Define the chemical potential and particle number as an expansion in A 



0- 



(2) , (3) , 

y^b = ^b + ^b + 



nb 



(2) , (3) , 

"-b + "b + 



(4.3) 
(4.4) 



The superscript denotes powers of Aq. In chapter 3 the results 



(2) 
^b 



n 



(2) 



aobmAQ 



am^ Aq 
Svr '' 



(4.5) 
(4.6) 



were found. Below, the next two orders in the chemical potential will be calculated. 
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It is useful to consider the molecular propagator at frequencies to <^ 1/ma? and 
momenta p <^ a^^. In the zero density approximation it becomes a free propagator 

Do{p, uj) ^ ^- — -, {p^a-\u;^ l/ma^) . (4.7) 

uj — p"^/2mb + ^0 ' 

In the first non-trivial order of the low density expansion, the normal and anomalous 
propagators were similarly seen in section 3.4 to be approximately Z times the usual di- 
lute Bose gas normal and anomalous propagators. The poles of the iterated propagators 
at low frequencies and momenta may then be separated as [38] 

u V 

Dn{p,uj) ~ Z T^V^^-^ ^ f\ 7T (p«a-\t^«l/ma2),(4.8) 

UJ — €[p) + lU UJ + e[p) — lU 

D.{p,uj) « Z^^^ + Z—^f^—, {p«a'\uj« 1/ma'), (4.9) 

UJ — e[p) + lO UJ + e[p) — lO ^ ' 



with 



^ 2e{p) ^ ' 



Here 



is the spectrum of low-lying excitations found by solving the Hugenholtz-Pines relation 
(3.39) at low frequencies and momenta. 

For the purpose of computing the vacuum scattering lengths a^f a-^d Ob in chapter 
3 there was only one relevant length scale, namely the atom-atom scattering length a. 
However, in the BEC regime another length scale is given by the average interparticle 
spacing n''^'^ and na^ is a dimensionless parameter. It is now seen how another impor- 
tant length scale appears of order -j^ ^ a. This scale, known as the healing length, 
defines a crossover in the excitation spectrum of low-lying excitations; for p ^ ^/na 
these behave as sound waves with a constant velocity while for ^/na <^ p <^ a~^ they 
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are almost free particles. Indeed, for small momenta 



/ (2) 

e(p)~P\/ — , p<.^/^ (4.13) 

V "ib 

as found in Eq. (3.55), while for higher momenta 

*>"2l + ''^"-(''i.")'^' v/iH«,«a-'. (4.14) 

Thus there are three regimes of momenta relevant for the problem, small momenta 
p <^ \fna^ intermediate momenta ^Jna <^ p <^ a^^ , and large momenta p S> a^^. In 
the first two 



li? ~ f ^ ~ ^ ) P <. y/na, (4.15) 

i" f 2p 

n;«l, t-p = ^ — 4^, ^Jna<.p<.a~^. (4.16) 

In the limit of large frequencies and momenta, the square root structure of the molecular 
propagator dominates, and in this limit it will no longer be possible to separate the poles. 

4.2 Corrections to the chemical potential and ground state energy 



In this section the first corrections to the chemical potential and ground state 
energy will be calculated in the low density expansion. This will be done by a very direct 
method, utilizing the particle number equation for the fermionic atoms, Eq. (3.46), and 
the Hugenholtz-Pines relation, Eq. (3.44). 

The chemical potential of the bosonic molecules was found in Eq. (3.44) by 
requiring a gap-less sound mode in the BEC. The result was 

/ib = ^ [Sn(0) - S2o(0)] - ^^^ [Sn(0) - S2o(0)]' . (4.17) 

The lowest order contribution to ^b is of order Aq. The two non-trivial corrections to 
the chemical potential considered in this chapter are of order Aq and Aglog Aq. Thus 
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(a) 



(b) 



Figure 4.2: The normal self energy at order Aq. In this chapter, the thin wavy lines 
are molecular propagators Dq, while the thick wavy lines are the iterated molecular 
propagators D^ and Da [Eqs. (3.36) and (3.37)], distinguished by the directions of 
arrows. 



it is seen that the last term in Eq. (4.17) does not influence these corrections to the 
chemical potential and the Hugenholtz-Pines relation becomes 



^b~^[Sii(0)-S2o(0)]. 



(4.18) 



The difference between the normal and anomalous self energy of the bosons is 
needed at order Aq . The diagrams contributing to the normal self energy in this order 

(2) 

are depicted in Fig. 4.2. The diagram in Fig. 4.2a has no analogue in S]^^ . Its value is 



2i lim 



duj 



^trjuj 



cfip 



t^^Z)n(p,t^), 



(4.19) 



n^o+J 2tt~ J (27r)3 

where the factor 2 appears as a combinatorial factor and the exponential factor ensures 
the correct order of field operators. The four-momentum dependence of t has not been 
made explicit. The t-matrix has the property that it is roughly constant at oj ^ l/ma^ 
and p <^ a^^, while it drops quickly to zero for large frequencies and momenta. The 
value of the integral may then be estimated by assuming that the main contribution is 
from frequencies and momenta small compared with \Ei,\ and a^^. Using Eq. (4.8) for 
the normal propagator, the frequency integration is performed by closing the contour in 
the upper half plane. The integration picks out the part of the propagator proportional 
to Vp which is suppressed at momenta p ^ \/na, and the value of the diagram is 



x^ 






Figure 4.3: The anomalous self energy for the bosons at order Aq. 
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approximately 



+bb 



d'p -e(P) + 2k+^b 



(2) 



(4.20) 



(27r)3 e{p) 

Using the asymptotic forms (4.13) and (4.14) of the excitation spectrum e(p), the inte- 
grand is seen to go as p for p <^ ^Jna and as p~^ for ^/na <C a^^ . This means that the 
integral is dominated by momenta of order ^Jna -^ oT^ and the assumption that t 
could be taken constant is justified. The diagram in Fig. 4.2a then has the value 



16 iaab"iAQ /a^ am^ Aq 



vr Svr 



(4.21) 



(2)u 



The diagram shown in Fig. 4.2b is already contained in S]^^ if the iterated normal 
propagator is replaced by Dq- To compute the contribution to ^^l from the diagram, 
the part proportional to Aq must be subtracted and the value of the diagram is 



4iA2 



duo 



(fp 

2tt J (27r)3 



t^^ [D^{p,u;) - Do{p,u;)] Do{-p, -cv) 



16Z 



_iaah'mAQ /aj^ am^ Aq 



4 y vr Svr 

Again the integration is dominated by momenta of order p 
normal self energies gives 



(4.22) 
na. The sum of the 

ZE<?(0)=(f + 16)^^i2!^. (4.23) 

(3) 

The only contribution to the anomalous self energy T,\-^ in this order is given by 
the diagram shown in Fig. 4.3. 



duj f d^p 



2-K J (27r) 



t'^^Z)a(p,w). 



(4.24) 




->■ 
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Figure 4.4: A contribution to the anomalous self energy already included at order Aq. 



Assuming that the integration is dominated by small p, lo results in 

d^p 1 



Ztbb^(2) 



(27r)3e(p)' 



(4.25) 



which is formally divergent at large p. However, the divergence is proportional to Aq 

(2) 

and arises because the diagram shown in Fig. 4.4, and belonging to S-[^ , is included 
in the diagram in Fig. 4.3. The diagram should be subtracted and the anomalous self 
energy at order Aq becomes 



zs(J(o) 



Z2tbb^(2) 



(fp 

(2^ 



1 



e{p) 



2mb 
p2 



^aahmAl /a^ am^ Aq 



4 ^ IT 8tt 
With the regularization, the diagram is dominated by p ~ ^/na. 
The chemical potential up to order Aq is found to be 



(4.26) 



Atb 



Z[Sn(0)-S2o(0)] 



aahmAQ 



1 + 



40 alam'^A?, 

/ V Svr 



(4.27) 



(3) 

The particle number must be determined in the same order, n^ originates from 
only one diagram, shown in Fig. 4.5, with the value 



n 



(3) 



1 



-2i^ lim 



v~ 



dpo 



0+ / 2-K 



-2Z lim 



»;^0+ J 27r 



dpo 



flWO 



,ITIP0 



^'^ ^ ^'^ [G{p)]'G{-p + q)D^{q) 



(2^4) 



1 



d\ 



(2^)3 

d^p 

(27r)3 \pQ - p^/2m - l/2ma'^ + iO^ 

1 qo + gV2mb + /u[,^^ 



(27r4) -po + go - (-P + g)V2m - l/2ma2 + iO 



e(g)2 + i0 



(4.28) 
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Figure 4.5: The diagram giving the particle number at order Aq. 



where it has again been assumed that the integration over q is dominated hy q <^ a^^ . 
The factor ^ comes from nb = f^/2 while the factor 2 is from the sum on fermion spins. 
Integrating over frequencies and shifting the momentum p -^ p+q/l to get rid of angular 
dependence, the expression reduces to 

^(3) ^ z f ^^( ^ y -e(g[+gV2mb+^ 

^"^ J {2'nY{2TTf\e{q)+p^/m + q^/2m^ + l/ma?) e{q) 

The q integration is indeed dominated by momenta q of order ^/na. The integration 
also results in terms of order Aq and higher which have been ignored. 
The particle number to order Ag is thus found to be 



Combining this result with the result for the chemical potential, Eq. (4.27), results in 



.. = ^^^^ 1 + fJ^ , (4.31) 

"lb I 3 V TT / 

exactly coinciding with the result in the standard dilute Bose gas [45]. 
The chemical potential is related to the ground state energy by 

I dEo , , 

"" = vai- *"■''' 

Assuming a form Eq = an^ + Pn-^ the ground state energy is found to be 



^0 27Taknl 128 /ggnb ■ 
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As this result is computed from the chemical potential it of course also matches the 
usual dilute Bose gas result. 

4.3 Next order corrections 

An elegant derivation of the chemical potential and ground state energy in the 
standard dilute Bose gas was presented by Hugenholtz and Pines in Ref. [38]. These 
authors noticed that the two quantities are linked by two equations. The first of these 
was given as Eq. (3.49) and is 



''"'"^ = h ,5SV / ^''"^ / W^ (^ +^V2mb) A.(P,^), (4.34) 



while the second equation is Eq. (4.32). The crucial observation was that given the 
normal bosonic propagator in some order in the low density expansion, when inserted 
in Eq. (4.34) the resulting relation between the chemical potential and the ground state 
energy will be an equation in the next order. This will be demonstrated below. 

The method may be used to compute the results of section 4.2 above. Inserting 
the propagator D^ from Eq. (4.8), correct up to order Aq in Eq. (4.34), results in 

- - ^nbMb = -- lim / -e^"- I /^ (u^ + p'/^m^) f, , (4-35) 

V 2 ^^^ 2n^o+J 2tt J (27r)3 ^ ^' ^cj + e(p)-iO' ^ ' 

where it has been assumed that only small p, uj contribute to the integral. Carrying out 
the integration, the assumption is seen to be valid, and the result is 

V 2 lovr^mb V '^2tt 



2 



-{Ana^f/^nl/'. (4.36) 



ISvr^mb 

In the last step, Ag has been traded for nb using Eq. (3.47). This is correct at the 
current order in na'^. Now assume the following forms for the ground state energy and 
chemical potential 

7^ o 2 5/2 5/2 

En 2TTaunf a/ n. 

-f = ^^ + a^^ ^, (4.37 

V TTT-b f^h 
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. 5/2 3/2 

Mb = ^^^^^ + /^ "^ ■ (4.38) 

mb mb 

Performing the differentiation in Eq. (4.32) gives 

5/2 3/2 

Mb = ^^^^ + ^a^^^^^, (4.39) 

from which it is seen that /3 = |a. A second hnear equation is provided by comparing 



powers of n^ in Eq. (4.36). Solving the equation gives a = — r^ ) a-nd inserting 



a 



and /3 in the equations for ^b and i?o, the resulting ground state energy and chemical 
potential are seen to exactly match the results in Sec. 4.2 above. 

4.3.1 The logarithmic contributions 

In the low density expansion, the next order contribution (computed by Wu and 
Sawada [76, 65] and denoted below by WS) to the chemical potential and ground state 
energy is logarithmic, proportional to log -^^. This in turn means that the correction 
arises from momenta in the range ^/rla <^ p <^ a^^ . It should thus not come as a great 
surprise that the BEC of weakly bound molecules also at this order coincides with the 
standard dilute Bose gas since the precise structure of t will not be probed by these 
momenta. 

As above, the normal molecular propagator is needed. Define 

D^{p,Lo) = D^^\p,u) + D^;^\p,u) + ... (4.40) 

The propagator Da {p,u>) is the one considered above in calculating the Lee-Huang- 
Yang (LHY) ground state energy (4.33), and is correct to order Aq. D^ {p,uj) is the 
part of the molecular propagator leading to the logarithmic corrections while the dots 
correspond to even higher order parts of the normal propagator. The anomalous prop- 

f2) 

agator Da {p,^) is defined in the same manner. 

In Eqs. (3.36) and (3.37) the normal and anomalous molecular propagators were 
related to the "bare" propagator Do and the exact self energies, the result of Beliaev 
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Figure 4.6: The expansion of the normal propagator, Eq. (4.43). The circles are the 
effective self energies Ti'^^ and S2o- 



[8]. It is possible to derive an analogous expression relating the exact propagators not 
to the bare propagator but to Dn {p, (^) and Da {p, ^) [38]. This requires the definition 
of new effective self energies 



^[^{p,L0) = ^ll{p,Lo)-^'SiO,0) 



(4.41) 
(4.42) 



(2) (2) 

since the self energies (at zero four-momentum) 'E\^ and S20 are already included in 
the LHY propagators D^ and Da • In order to compute the WS order, it is sufficient 
to approximate the normal molecular propagator by [38] 

D(2) (p^ ^) = Dip (p^ ^)s;^ (p, u;)D(^) (p, a;) + D,« (p, c^)S'2o(p, u;)dP {p, u) 

+d(^\p,co)^',o{p,u;)d(^\p,lo) + D(^\p,Lo)^[,{-p,-u;)Di'\p,u;). (4.43) 

Diagrammatically, the expansion is shown in Fig. 4.6. Because of the range of momenta 
relevant to the problem (p ^ ^/na) , the propagators may be expanded in powers of the 
density and become 



dW(p,u) « Do(p,^) « 27o T^' 

uj — p"^/2mb + ^0 

dI^\p,uj) « Alt'^'^Do{p,Lo)Do{-p,-uj) « 



^/^f 



to — p'^/2mb + iO —UJ — p^ /2m\y + iO 

(4.44) 

(2) 



Inserting the resulting normal molecular propagator D\i in Eq. (4.34) results in 
an equation between the ground state energy and the chemical potential in WS order 



Eo 1 
V 



iZ 



-nhl^h 



lim 



duj ,-, 



trjui 



0+ ./ 2-K 



d?p 

(27r)3 



i:ii(p,u;)((j + p^/2mb) 
UJ — p^ /2m\y + iO 
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Figure 4.7: The class of diagrams which must be considered in the anomalous self energy. 



(2) 






+ 



^'ii{-P,-^)fJ'h^ 



. (4.45) 



The quantity in square brackets is {lo + p^ /2m\^)D\i {p,u>) and the normal propagator 
has the property that any poles must lie slightly above the negative real uj axis or slightly 
below the positive real lo axis. Looking at the a; integration, which must be closed in 
the upper half plane, it is then concluded that the following contributions to the self 
energies need to be considered: 



• poles of S2o(^''^) above the negative real axis 

• poles of Ti'-^-^{p,uj) above the negative real axis, and 

• poles of Ti'-yT^{—p^—Lo) above the negative real axis as well as the momentum 
dependence of S^-^ (p,p^/2). 

The simplest class of diagrams entering the anomalous self energy sum to give the 
molecule-molecule scattering t-matrix. This reduces to a constant at p ^ a~^ and this 
class of diagrams does not have any poles at relevant momenta. Poles above the negative 
real exis become possible if, in the molecule-molecule scattering problem, a pair of bare 
molecular propagators are replaced by the iterated propagators D^^\ thus all diagrams 
of the form depicted in Fig. 4.7 must be taken into account. The leading of these 
diagrams is shown in Fig. 4.8. Carrying out the frequency integration, anticipating 
that the main contribution will be from p <^ a^^ and uj <^ 1/ma^, results in 

2 f d'^q rrib 



E'2o(p,^)=4A2zVr * 



2 -72,, (2) f.hh 



(27r)3 g2 (^ + g2/2r7ib + /c2/2mb - iO) ' 



(4.46) 
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Figure 4.8: The diagram belonging to Sgolp, tj) which gives the leading contribution in 
Eq. (4.45). 



For simplicity k = p + q. The tilde signifies that this is the only relevant part of 
the self energy at WS order. In the prefactor, the factor 4 comes from combinatorics, 
Z'^fi'ly from propagators, Aq from condensate lines, and [t ) from the two t-matrices 
evaluated at low momenta. 

The normal self energy is constructed in the same way. The two leading diagrams 
relevant for computing the poles of T,ii(p,uj) above the negative real uj axis are shown 
in Fig. 4.9. The sum of these diagrams is 



t[,ip,.) = -4AlZ^U^^Y(t-^'' 



1 1 

+ 



mt 



d?q 
(27r)3 V^g2/.2 ' /c4y ^_^g2/2mb + A;V2mb-iO' 

(4.47) 



where it is again anticipated that only momenta small compared with the inverse scat- 
tering length are relevant. 

The final contribution from molecular self energies comes from S']^^(— p, — w) and 
here it is important to include the momentum dependence of the self energy. This 
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Figure 4.9: The diagrams belonging to ^'-^^(p, cj) which result in the leading terms in 
the low density expansion of Eq. (4.45). 
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Figure 4.10: A contribution from T,'^^(—p, —uj). 



may be included by considering the relation between the normal self energy and the 
scattering amplitude. This relation was first derived by Beliaev [8, 9] (and similar 
results for fermion scattering by Galitskii [31]) and at relevant momenta it takes the 
approximate form [38] 



Sii(p,w) 



+2Z^Al (t'^b 



2Allmfsip/2,p/2) 
2 f (fa 



1 



+ 



.(4.48) 



(27r)3 \uj — p^ /Am\y — q^ /m\^ + iO cp' jm^, — p^ /Am\^ — iO^ 

Here fsip,Q) = ifip,Q) + fi~Q:P)) f^ is a symmetrized scattering amplitude and re- 
duces at small momenta to t . Note how the first term under the integral simply 
corresponds to the diagram of Fig. 4.10 while the second term under the integral is a 
kinematic term. The reasoning behind Eq. (4.48) is that while the scattering amplitude 
between two particles contains the possibility of multiple scattering, it is still an intrin- 
sic two-body quantity in the sense that it is assumed that after the scattering the two 
particles are in the vacuum and on-shell which is certainly not true in the many-body 
case. The proof by Beliaev of this equation is for a dilute gas of bosons but it goes 
through in exactly the same way in the problem of the weakly bound bosonic molecules, 
since the assumption is that the potential between the bosons (which in the case of 
the fermion problem is replaced by the two-boson irreducible diagram F, see chapter 3) 
reduces to a constant at momenta p <^ a^^. 

The imaginary part in the first term of Eq. (4.48) may be written in a different 
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7^ -^s 

(a) (b) 

Figure 4.11: An illustration of how fermions may be thought of as constants at low 
momenta. Fermions lines are drawn as straight thin lines. 



way, as also found by Beliaev, 

Im/,(p/2,p/2) = -Z^ (t'^b)'lm / yf^^^ \^, -■ (4.49) 

The real part of the first term in Eq. (4.48) is simply the scattering matrix t which 
reduces to a constant and does not contribute to the ground state energy. 

At this point it is important to investigate whether the presence of the fermionic 
atoms affects self energies at the current order. The effect of replacing a bare molecular 
propagator insde the two-boson irreducible diagram by an iterated propagator Dn 
must be examined. As an example, consider the diagrams in Fig. 4.11. The fermionic 
propagators in this diagram all have chemical potential n = —l/2mo? in the BEC regime 
of the gas, see Eq. (3.40). Equivalently, when a molecule in the vacuum breaks into two 
atoms they equally share the binding energy —1/mo?'. For momenta small compared 
with the inverse scattering length, the fermionic propagators are thus approximately 
constant and the diagram shown in Fig. 4.11a effectively shrinks to that of Fig. 4.11b. 
The integration over the normal propagator closed on itself was performed in Eq. (4.21) 
and goes as Ag. Thus the diagram in Fig. 4.11a goes as Ag and is of higher order than 
the WS order. The same argument goes through for all two-boson irreducible diagrams 
with molecular propagators replaced by iterated propagators. 

The fermionic atoms are also modified by their self energies. The anomalous 
fermion self energy at lowest order is simply two fermions coming in and a condensate 
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line going out. Such anomalous self energies must come in pairs and thus modifying 
the two-boson irreducible diagram by the simplest two such anomalous self energies 
simply corresponds to adding the possibility of three boson scattering. This in turn 
only contributes at order Aq. The normal fermion self energy consists of atom-molecule 
scattering. Replacing a fermion propagator with the atom-molecule scattering process 
inside the two-boson irreducible diagram again results in a diagram belonging to the 
class of three boson scattering. Furthermore, the atom-molecule diagram reduces to a 
constant at p <^ a^^. Thus the presence of this self energy does not change the pole 
structure of the two-boson irreducible diagram in the range of momenta and energies 
under consideration. 

Each two-boson irreducible diagram will contain exactly one fermion loop. The 
external four-momentum going into the two-boson irreducible diagram may be contained 
completely inside the loop and since the fermionic propagators are constant at p <^ 
a^^, the pole structure of the molecular self energies at small momenta is completely 
unchanged by the presence of fermions in the problem. 

The conclusion to this discussion is that the presence of fermionic atoms does not 
change the chemical potential and ground state energy from the results of the standard 
dilute Bose gas beyond replacing the scattering length Ob with 0.60a. 

The rest is mathematics. Inserting the self energies, Eqs. (4.46), (4.47), and 
(4.48) in Eq. (4.45) and performing the integrals, the result is [38] 



The range of integration is determined by p being in the range ^/na <^ p <^ a~^ and 
with logarithmic accuracy ^/rla and a^^ may be used as the actual limits of integration. 
Using the relation (4.32) between the ground state energy and the chemical po- 
tential as well as the particle number equation (4.30) the result of the usual dilute Bose 



gas [76, 65] is reproduced 

^ = levr^^LJl ^ ^ log nbO^ , 4.51 

H^ = ASn'^^l^-^^login^al). (4.52) 

It has been assumed that rib does not contain terms of order AglogAo. Arguments 
similar to those given above show that this assumption is indeed correct. 

In conclusion, the LHY and WS terms in the low density expansion arise from 
non-trivial renormalizations due to the infrared divergences of molecular propagators. 
At low momenta the fermionic atoms in the BEC regime are nearly constants and their 
dynamics do not affect the low-density expansion of the chemical potential and ground 
state energy at the order considered here. However, from the above discussion, and also 
from the Hugenholtz-Pines relation (3.44), it is obvious that the presence of fermions 
will affect next order quantities in a non-trivial way. 



Chapter 5 
Two-component Fermi gases with a mass imbalance 



Since the experimental discovery of the BCS-BEC crossover in ultracold mixtures 
of atoms in two different hyperfine states of e.g. ^"^K [69], a major effort has gone into 
creating other more exotic superfluids. An example is a system with a mass imbalance, 
where instead of studying two different hyperfine states of the same fermionic atom, a 
mixture of two different atomic species is used. While experiments in these mixtures 
are still in their infancy, these systems have received considerable theoretical attention 
[56, 59, 39]. The interest in these systems has several reasons. For instance, the advent of 
the use of optical lattices in ultracold gases has opened up the possibility of continuously 
changing the effective mass of one of the species independently of the other, since the 
resonant frequencies of the atomic species will be different. 

One might naively expect that the presence of a mass imbalance would serve only 
to modify the few-body coupling constants computed in chapter 3. However, it turns 
out that beyond a certain mass ratio, the physical picture becomes quite different. It 
was discovered in Ref. [57] that, as the mass imbalance is increased, collisional losses of 
the Fermi gas also increase. For small mass imbalances, the gas becomes more stable 
as a Feshbach resonance is approached. Conversely, beyond a mass ratio of 12.33 the 
system becomes less stable as a Feshbach resonance is approached and finally, beyond 
the mass ratio 13.6, the three body system was found to display bound states in a 
manner very similar to the Efimov phenomenon in systems of three bosons [25]. 
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For moderate mass ratios, the main effect of introducing a mass imbalance is 
indeed to change the few-body couphng constants. Results for the three and four fermion 
problems, calculated in coordinate space using the Bethe-Peirls method, have been 
obtained in Refs. [56, 58]. The results for four fermions have been confirmed by Monte 
Carlo methods in Ref. [75]. Thus, the results in this chapter for three- and four-body 
scattering calculated using the momentum space formalism are not new, rather they 
confirm already known results by the use of a complimentary method. 

Even though the Fermi energies will be different, in heteronuclear mixtures with 
equal densities of the two atomic species, the Fermi momenta will be equal. This means 
that the usual Cooper pairing is still taking place on the BCS side and the picture of a 
BCS-BEC crossover is valid. In this chapter the BEC regime with equal densities will 
be examined. 

5.1 Molecular propagator 

The formalism presented here is very similar to that presented in chapter 3. 
Consider two species of fermionic atoms interacting close to a heteronuclear Feshbach 
resonance. For simplicity of notation, even though the particles are not different hyper- 
fine states of the same isotope, the particles will be called j and |. Let the particles 
have masses 

m^ = 7m, nil = m (5-1) 

and chemical potentials ^|,/ij < 0. 

To calculate the molecular propagator, the polarization bubble which appears in 
the repeated scattering between the two species of atoms is needed (see Fig. 3.1). The 
polarization bubble has the value 



d?q 1 



(27r)3 LO - (g +p/2)2/2m| -{q- p/2)2/2m| + ^| + /i^ + iO 
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= 2~ + ^^^V~^+^ /2"^b-^T -A^i -«0, (5.2) 

where it has been used that q can be shifted such that the angular integration becomes 
trivial (let q^ q+ | ^|~^| p)- Here 

rrir = — ^— (5.3) 

m| + mi 

is the reduced mass and m\j = 171^ + ttij^ is the mass of the molecule. 

The scattering length between the two species of atoms is related to the vacuum 
polarization bubble (for which ij,^ = fi^ = 0) by 

rrir-^, . nir- —A / 2-k 2A\~ , ,, 

«=2^^<°' = 2^rTAn(0) = ("i^ + v) • <"' 

The molecular propagator is then given by 

D{p,uj) = ^, ^ = — ^ (5.5) 

1 + XU{p, Lo) mr a-i - ^2m^^J-uj + p^ /2m), - fi^ - m - iO 

which reduces to the usual form (3.16) in the limit of equal masses. According to 

Hugenholtz and Pines [38], the molecular propagator should have a pole at p,uj = 0. 

This leads to the relation 

^b = MT+/^i = -^ (5.6) 

for the binding energy of the molecule. 

In the limit of small momenta and energies, p <^ a^^ and lo <^ l/2my-a'^, the 

underlying structure of the molecule is not probed. In this limit the propagator reduces 

to 

D(p,uj)^ Ti • (5.7) 

^■^' ' uj-p^/2mi, + i0 ^ ' 

Z = —^ is the residue and the remainder is the free propagator of a particle with mass 

mb = r7T,| + m^i as expected. 
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Figure 5.1: The integral equation satisfied by the atom-molecule scattering t-matrix. 



5.2 



Three-body problem 



Assume an atom of type | scatters off the bound molecular state. The general 
integral equation for the atom-molecule scattering matrix with kinematics as shown in 
Fig. 5.1 is given by 



t{p,Po;q,qo) = -Gi{p+q+E)-i 



d*Q 

(27r)4 



G^{Q + q+E)G^{-Q)D{Q + E)t{p,po;Q,QG 



(5.8) 



E is here used both as the total energy and as the total four-momentum (0,E). The 
symmetry factors are the same as in the mass-balanced case. 

As in chapter 3, in the notation for the t-matrix the dependence on the total 
energy has been suppressed. Also, E < E\^ in order for processes in which the molecule 
and atom scatters into a final state of three free atoms to be forbidden. 

The frequency Qo can be integrated out by closing the contour in the upper half 
plane, where the only pole is the pole of G|(— Q). This results in 



t{p,po;q,qo) = -Gi{p + q + E) 



where Qq = —Q'^/2m^. 



d?Q 

(27r)3 



Gi{Q + q + E)D{Q + E)t{p,po;Q,Qo), (5.9) 



Consider first the s-wave scattering length which is proportional to t(d, 0;d, 0) 
evaluated at E = Ey^. The incoming momentum and frequency may be taken to vanish 
but the outgoing should be kept finite in order to have a solvable integral equation. The 
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result is the equation 

1 +|)^ 

2ttp Jo ^1 - Vl + g2 [1 _ ^2/(1 + ^)2] ''^l+p2+q2_ 2qp/ (1 + ^-l) ' 

(5.10) 

where momenta are measure in units of inverse scattering length. Here the integral 
equation for t(p) = t{0,0;p,—p'^/2m^) is solved with pQ = —p^/2m,i^ to make the de- 
pendence on momentum the same on both sides of the equation. This point has been 
discussed in detail in section 3.5. 

Having solved Eq. (5.10) at a certain mass ratio, to calculate the scattering 
amplitude each external molecular propagator has to be renormalized by the square 
root of the residue of the pole of the molecular propagator. The scattering amplitude 
at vanishing momentum is thus 

T(0) = Zt(0). (5.11) 

The relationship between the scattering amplitude and the scattering length is 

T(0) = — Obf, (5.12) 

with 77231. = ^ ^ ^ being the three-body reduced mass. The three-body scattering 



length is found to be 



"^3r t(0) 
Obf = ^a. (5.13) 



The result is shown in Fig. 5.2 and exactly matches the result of Petrov [56]. 

5.2.1 p-^vave scattering and collisional relaxation 

Up to now the only bound molecular two-atom state considered has been the 
weakly bound and large state of size comparable with the scattering length a. This 
state is the state appearing at a Feshbach resonance and is the highest ro-vibrational 
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Figure 5.2: Obf/a as a function of 7 = m-^/m 



T/'"-!- 



state. Additionally, in the two body problem there are tightly bound states of size i?e, 
the scale at which short distance physics becomes important. Through collisions the 
weakly bound molecules can fall into these deeply lying states and in the process a 
binding energy of order l/(mi?g) is converted into kinetic energy. The released energy 
is quite large and the colliding atoms will escape the system. This process is called 
collisional relaxation and is the process which determines the lifetime of the Bose gas 
of weakly bound molecules [57]. 

The rate at which the collisional relaxation occurs is conventionally written as 
either 



nb = -Orcinb 



(5.14) 



or 



nb 



-/^rell-b- 



(5.15) 



Both versions are natural; for the first, two molecules need to approach each other closely 
for the relaxation process to occur and Oj-d is the probability for this to occur. For the 
second I5^^-^ is the characteristic relaxation time. Below, the result will be presented in 
the form (5.14). 
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For equal masses of the constituent fermions the relaxation rate Orel decreases 
with increasing a as [57] 

Qrel OC — ^ ( — ^ ) . (5.16) 

m \ a J 
This ensures the stability of the Fermi gas in the vicinity of a Feshbach resonance as 

opposed to the case of Bose gases which become unstable close to a Feshbach resonance. 

The reason for the stability of the Fermi gas is that the relaxation process requires three 

particles to approach each other within a distance of order R^, while the remaining 

particle in the molecule-molecule scattering will typically be at a distance of order a. 

The probability for three fermions to approach each other is suppressed because of the 

Pauli principle, since two of the fermions will necessarily be identical. 

It was shown in Ref. [56] that the leading relaxation channel at large a is scattering 

through the p-wave channel. In general the relaxation rate may be written as 

Orel oc — ^ ( ^ ) . (5.17) 

m \ a J 

The power s is a function of the mass imbalance, s is related to the power-law behavior 

i/j ~ p"^^^ satisfied by (the p-wave part of) the three-body wave function in the range 

Re <^ p -^ a hy s = 21/ + 1. This can be seen from estimating the probability for the 

three particles to be within a distance of Rg of each other. This probability is 

j^d^pmpw \a) ■ ^"-^'^ 

The relation between the outgoing part of the wavefunction and (the p-wave part of) 
the t-matrix is 

/d J) f do /-* \ 

(27l¥ / T^^TAi^^P^ ^"^^ \^'^'^P'P0) 

xG^{p -q + E)G^{-p)Gi{q)e''''<P^+P^'>+'''^-P^-'^'^-P\ (5.19) 

This relation is illustrated in Fig. 5.3. If the t-matrix scales as t{p) oc p at large 
momenta, p ^ o,^^, then 

s = -3 - 25. (5.20) 
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Figure 5.3: Relation between the scattering matrix and the scattered wave function. 



The colhsional relaxation rate may be calculated by noting that the scattering 
t-matrix can be written as the sum over partial waves 



tip,Po;q,qo) = ^(2^ + l)P£(cos0)tW(p,po;'7,' 



(5.21) 



£=0 



Here 6 is the angle between p and q. Pi (cos 9) are the Legendre polynomials of which 
the first few are 



Pi(x) 



1, 



X, 



P2{X) = ^{ZX^-I). 



These satisfy the orthogonality relation 



1 2 

Pi>{x)Pi{x)dx= ^j^' 



(5.22) 



To see that the integral equations for the partial waves decouple, consider the 
general integral equation, Eq. (5.8). Let 9 denote the angle between p and q and 
likewise (p the angle between q and Q. Then the angle 7 between p and Q satisfies 
cos 7 = cos^cos(/) + sin ^ sin cos A. With this definition of angles, only the t-matrix 
under the integral depends on the azimuthal angle and the corresponding integration 
may be carried out as 

f.27r 
/O 



/ —t{p,po;Q,Qo) = y2PE{cos9)Pe{cos^)t^'\p,po;Q,Qo). (5.23) 
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Here the addition theorem for Legendre polynomials [44] 

Pi{cos-f) = Pe{cos9)Pe{cos(j)) + ^ 2|- ^-Pp{cose)Pp{cos(l))cosmX (5.24) 

m=i 1^ + "^)- 

has been used. The integral equations are now decoupled by using the orthogonality 

relation 5.22 with the result 

1 f^ 

"i^ / Q^dQt^^\p,po;Q,Qo)D{Q + E) d(cos</.) PKcos(/>)G|(g + Q + ^), (5.25) 

where Qo = — (5^/2m|. 

To solve the integral equation for the p-wave contribution, let qq -^ —q'^/2mp 
For simplicity let also pq -^ —p^/2m,i^. The p-wave contribution to the t-matrix satisfies 

t'P{p, -p'^/2mf, q, -q^/2m^) 

2(1 (l+7~^)(l+P^ + g^) i l+p^ + q^+2pq/{l+^-^] 
— Tfid 1 /I 9 9 ^^-'S 



pq Ap^q^ 1 -\- p'^ -\- q^ — 2pq/{l +7 ^] 

, 1 + 7'^ r Q2^Q_t{P^-pV'^rnY,Q,-Q^/2m^) 



^ Jo "1- VI + Q2 [1-^2/(1+^)2] 

/^ _ (l+^-l)(l + g2 + Q2) l+g2 + Q2+2gQ/(i+^-l) X 

VgQ 4g2Q2 °g 1 + ^2 + g2 _ 2<^Q/(i + ^-1) ; ^^- ""> 

at E = E\^. Fig. 5.4 shows the dependence of s on the mass imbalance found by solving 
Eq. (5.26). The curve matches the result of Ref. [58]. 

It is seen that above a mass ratio of m-^/ni]^ ~ 12.3 the relaxation rate increases 
with increasing scattering length. This is nicely explained in the Born-Oppenheimer 
picture [58]; in the case of a large mass ratio, if the two heavy atoms are separated 
by a distance much smaller than a, then the presence of the light atom will mediate 
an effective attractive -^ potential between the heavy atoms. On the other hand, the 
Pauli principle results in the centrifugal barrier, which is a repulsive ^ potential. The 
mass imbalance of 12.3 corresponds to the cancellation of these potentials in the Born- 
Oppenheimer picture. As the mass ratio is further increased, the attractive potential is 



s 1 




Figure 5.4: The power s appearing in the relaxation rate, Eq. (5.17), as a function of 

7- 



also increased and at a mass ratio of approximately 13.6 one encounters the phenomenon 
of the fall to the centre [44] : The system displays a behavior characteristic of the Efimov 
phenomenon [25] where a hierarchy of three-body bound states appear [58]. The short 
distance physics becomes very important and a short range three-body parameter is 
needed in order to accurately describe the physics, see e.g. Refs. [7, 12]. 

5.3 Scattering of composite molecules 



Consider now the 4-body problem with a mass imbalance. The present calculation 
follows closely the mass-balanced molecule-molecule scattering calculation presented in 
section 3.6. Thus some intermediate steps in the calculation will be skipped. 

The governing integral equation is formally identical to the problem of equal 
masses, Eq. (3.67), that is 

i(p,Po)=r(0,0;p,po) + ^ / q'^dqdqoD{q,qo+E/2)D{q,-qo + E/2)r{q,qo;p,po)t{q,qo). 

(5.27) 
The difference from the previous calculation is that now both the molecular propagator 
and the two-boson irreducible vertex, T, depend on the mass-ratio. 
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Below, the convention is used that all atomic propagators have their frequency 
shifted by a quarter of the energy going into the scattering, while the molecules have 
their frequency shifted by half the energy. That is 

Po — p /2m'i^l + E/A + iO 

D{p,po) = — ^ (5.29) 

rrir a^i - ./2rnry/-po - E/2 + p2/2mb - iO 

The relation between T and H is 
T{q,qo;p,po) = \j (^G't(p/2 + Q)G^{p/2 - Q)E{q;p/2 + Q;p/2 - Q), (5.30) 
and the vertex H satisfies 

E{q;p/2 + Q;p/2-Q) =- j '^ [G^{-p/2 - q + Q)G^{-p/2 + q-Q) + {q^ -q)] 

-i j ^ {Gi{p/2 - Q')Gt(-3p/2 + Q')D{-p -Q + Q')E{q-p/2 + Q-p/2 - Q') 
+Gt(p/2 + Q')Gi(-3p/2 - Q')D{-p + Q- Q')E{q;p/2 + Q';p/2 - Q)} . (5.31) 

The vertex E{q;p/2+Q;p/2—Q) may again be split into parts E'^{q;p/2+Q;p/2 — 
Q) \E^ {q;p/2 + Q;p/2 — Q)] analytic in the upper [lower] half planes of Qo, with 

S(g;p/2 + Q-p/2 - Q) = r + {q;p/2 + Q-p/2 - Q) + E-{q-,p/2 + Q-,p/2 - Q). (5.32) 

To achieve this use Eq. (3.73). Now insert Eq. (5.32) into Eq. (5.30) and perform the 
frequency integrations over simple poles (fermionic propagators). The iterated terms 
in the integral equations for H^ evaluated at the values of frequency dictated by the 
simple poles will then only depend on an on-shell vertex 

S(9;Pi,Pi/2rnt - E/4;p2,pl/2mi - E/A) = x{q;Pi,P2)- (5.33) 

Again, as in section 3.6, the on-shell frequencies would have been obtained had G^{pi) 
and Gi{p2) been attached to the vertex H(g;pi;p2) and their frequencies been integrated 
out using a contour around their simple poles. 
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The solution of the integral equation for the on-shell vertex x is now possible with 
the substitution 

po -^ {p/2 + Qf/2m^ + (p/2 - Qf /2mi - E/2. (5.34) 

The equation for the on-shell vertex is 



f dilff 

x{q;Pi,P2) = - ^[Gi{q-pi)G'i{-q-p2) + {q^-q)] 

- / 7^;^ S Gii-pi -p2 - Q)D{-Q - pi)x{q;pi,Q) 



o=Q2/2m|-£;/4 



+ Gi{-p,-p2-Q)D{-Q-p2)x{q;Q,P2) ^ ^^^^ ,,JA5.S5) 

where (pi)o = Pi/2m^ — E/A and (^2)0 = pl/^^^i ~ E/A. 
The relation between T and x is 

1 [ d^pi d^p2 



2 I (27r)3(27r)3 t^i^^ ~ ^i)^t(-P " :P2) + (p ^ -p)\D{-pi - P2)x{q;Pi,P2)- {5M) 

Again (pi)o = Pi/2m-^ — E/A and (^2)0 = Pi/2wij. — E/A and po is a free parameter. 
r*-^' is the result of calculating the Born diagram, 

rW(g,go;p,Po) = -iJ^J^G^iQ+p/2 + q/2)G^iQ-p/2-q/2) 



X 



Gii-Q-p/2 + q/2)Gi{-Q+p/2 - q/2) 



dn^ f d^Q A 



'^ ' AtT J (27r)3(^2_52)(^2_C'2)' (^"^^^ 



with 



A = E/2 - Q^ /2mr - p^ /^m^ - q'^/SrUr - p- q/Am^ +p- q/Anii, (5.38) 
B = po-Q- p/2m^ + Q ■ p/2mi - Q ■ q/2m,, (5.39) 

C = qo-Q- q/2m^ + Q ■ q/2mi - Q ■ p/2m,. (5.40) 

Finally, to relate the scattering amplitude to the scattering t-matrix at vanishing 
incoming energy and momentum use 

r(0) = Z^t(O), (5.41) 
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Figure 5.5: The scattering length at as a function of mass-imbalance, 7 = m-^/mi. 
Above a mass ratio of about 7 ~ 10 convergence becomes very slow and numerical 
errors quite large. The curve is thus only plotted to 7 = 12. 



while the scattering amplitude is related to the scattering length by 

r(0) = -a. 



(5.42) 



Mr = -Y^TTi is the reduced molecular-molecular mass. The result is shown in Fig. 5.5 



and matches the result of Refs. [58, 75]. 



Chapter 6 
Superfluidity close to a p-wave resonance 



Attention will now be turned to atomic Fermi gases interacting close to a p-wave 
Feshbach resonance. At low energies, because of the angular momentum barrier, typi- 
cally s-wave scattering will dominate. However, in gases consisting of identical fermionic 
atoms, the Pauli exclusion principle prevents s-wave scattering. At low temperatures in 
these systems, p-wscve scattering will dominate over higher angular momentum scatter- 
ing, and the study of p-wave resonantly paired Fermi gases is therefore of fundamental 
interest. 

The p-wave resonantly coupled superfluid state has not yet been observed ex- 
perimentally in degenerate atomic gases, but p-wave Feshbach resonances have been 
observed in both ^"K [69] and ^Li [67]. Considerable theoretical attention has been 
given to the subject [36, 22, 35, 17, 34], particularly since it was pointed out [35] that 
the two-dimensional p-wave superfluids may display phases which are candidates for 
quantum computing systems. Below, only three-dimensional systems will be consid- 
ered. 

The gas of p-wave resonantly paired fermions is in many ways a richer system 
than the s-wave system considered in the previous chapters. In the s-wave paired 
superfluid, the BCS and EEC phases were merely different limits of the same phase, 
but this is no longer true in p-wave systems. The order parameters now correspond 
to different projections of angular momentum and distinct symmetries. This allows for 
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the possibility of genuine phase transitions which in some cases can even be topological 
[74, 62, 22, 35] and the phase diagram becomes much richer. 

An interesting theoretical feature of the p-wave Feshbach resonances is that they 
are naturally narrow, whereas most s-wave resonances studied in experiments are wide. 
This means that a quantitatively correct perturbative approach may be used across 
the whole crossover [34]. However, as discussed below, the Feshbach resonances fall 
into categories of weak and strong, and while the weak resonances are amenable to a 
mean field treatment, as performed in Ref. [34], the strong resonances must be treated 
more carefully. In particular, it will be shown in chapter 7 that the system close to a 
strong resonance will contain three-body bound states (trimers), an effect first noticed 
by Y. Castin and co-workers [16]. 

In this chapter, the two-channel model for p-wave Feshbach resonances will be 
introduced and its parameters related to those of basic scattering theory. Next the 
important distinctions between wide and narrow and between weak and strong Feshbach 
resonances will be discussed. 

6.1 A two-channel model for p-wave resonances 

Consider a system of N identical (atoms in the same hyperfine state) spin-less 
fermions of mass m. At low energies, the scattering will proceed entirely in the p- 
wave channel, since s-wave scattering is excluded due to the Pauli principle and higher 
order scattering is suppressed. As in the s-wave problem, it is convenient to describe 
the system by a two-channel model. For the p-wave resonantly coupled superfluid the 
two-channel Hamiltonian is given by [71, 70, 37, 35, 17, 34] 



(6.1) 
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The operators at, Op- are creation and annihilation operators of the fermionic atoms 






of momentum p. Similarly, b\^ ^, 6^ p> are creation and annihilation operators of the 



bosonic spin 1 molecules with the vector index fi being the spin projection on some 
axis. The coupling 5(|p|) appears as the amplitude for the transition from a pair of 
identical fermions with relative momentum p and orbital angular momentum 1 into the 
bare closed-channel molecular state with internal angular momentum 1. 

The superfluid described by the two-channel Hamiltonian (6.1) is controlled by 
four parameters. The first of these is the bare detuning eo. This may in principle have 
a directional dependence, however for simplicity it will be taken constant in this thesis. 
The second parameter is the particle number N, which is the expectation value of the 
operator 

^ = E4^P^+2E^lAr (6-2) 

It is convenient to trade the particle number for the energy scale 

ep = (^-^ J /2"^- (6.3) 

In the absence of the Feshbach molecules this would have been the Fermi energy. The 
remaining two parameters are contained in the coupling constant fi'dpl). Taking g con- 
stant would imply that the interactions take place at a point in space, as if the molecules 
were infinitesimally small. However, in practice the molecules have a finite size set by 
the range of the interactions, i?e- fi'di'l) is in turn proportional to the wave function 
of the molecule in momentum space and should be taken approximately constant for 
p ^ A ~ p- while quickly dropping to zero for p ^ A. Thus define a cut-off function 

9{m^9iM). (6.4) 

In this introductory chapter to the physics of p-wave resonantly coupled superfluids, the 
simplest cut-off function, a step function, will be used. This choice is 

C(|p1) = e(l-|p1/A), (6.5) 
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which is 1 for |p| < A and above. In chapter 7 both this and a "softer" cut-off wih be 
investigated. 

6.1.1 The diagrammatic approach 

The propagator of fermionic atoms is 

G(p,uj) = ^ . (6.6) 

From the Hamiltonian (6.1), the propagator of the bare closed-channel molecule may 
be read off to be 

Do^^u{p,^) = ^. — T7. (6.7) 

cj — eo — p'^ /Am + i\) 

As in chapter 3, the bare molecular propagator needs to be corrected by self energy 
insertions as shown in Fig. 3.2. The physical dressed propagator is then 

DApM={-^-r, ^, r) • (6.8) 

The value of the self energy bubble is 

T.^,i,{p,uj) = 

2 / dqQ d?q qfiQu^.'^M) 



2ig' 



2tt (27r)3 [uj/2 + go - {q + p/2)^/2m + iO] [lo/2 - Qo - (q - p/2y/2m + iO] ' 

(6.9) 



The factor 2 appears as the appropriate combinatorial factor; in the s-wave problem, the 
operators giving the propagators of the two distinguishable fermionic atoms inside the 
bubble could be contracted in only one way, in the present case there are two possible 
ways. Notice how the physical molecular propagator will in general depend on the 
chosen method of cut-off. To proceed, the cut-off function needs to be specified. Using 
the simplest cut-off, Eq. (6.5), the self energy becomes 

2 2 1/"^ «4 



3 2tt'^ Jq uj — q^ /m — p^ /Am + iO 



Sfj.u { -ci + C2(g'^/4m - qo - iO) 



-C2^{q^/4m -qo- zOf/^^n-^' ^ . (6.10) 

A vg^/4 - mqo - lO j 

The integral scales as g^ for large momenta. Thus, two cut-off dependent parameters 

have been defined, 

^1 = ^^' ^2 = ^^. (6.11) 

c\ has dimensions of energy while C2 is a very important dimensionless parameter as 
will become clear below. 

The propagator of bosons becomes 

D^,y{q,qo) = 



(1 + C2)(go - g2/4m -UJO + iO) + C2^{q^/4m - qo - zO)3/2tan- 



-1 A 



■^g2/4-mi3o 

(6.12) 



which is diagonal in spin indices. It is seen that ci disappears from the theory, it is 
absorbed in the physical detuning^ as 

c^o = ^^^^ < 0. (6.13) 

1 + C2 

UJo is an important parameter in the theory; a negative ujo implies the presence of a 
physical bound state in the two-body system. On going through ojo = the bound 
state disappears and is replaced by a resonance. 

6.1.2 p-wave scattering 

At low energies, the scattering amplitude between two identical fermionic atoms 
may be approximated by [see Eq. (2.4)] 



^ In the two-channel s-wave model, a coefficient proportional to the cut-off appeared to shift the bare 
detuning to the physical detuning, see Eq. (2.16). No analogue of C2 is present in the s-wave problem. 
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with k being the relative momentum of the two scattering atoms. Whereas the s-wave 
scattering amphtude fs{k) -^ —a as k ^ 0, fp{k) ^ in the low-energy limit, demon- 
strating how s-wave scattering, if allowed, would dominate the scattering problem. The 
scattering volume v is the analogue of the s-wave scattering length, diverging and chang- 
ing sign across the Feshbach resonance. The parameter k^ is a characteristic momentum 
and is analogous to the effective range tq. 

In the model above, with the choice of cut-off function Eq. (6.4), the scattering 
amplitude is [35, 34] 

fv{k) = — ) -T , (6.15) 

from which it is possible to identify 

mg'^ , 4A(1 + C2) 

67r(l -I- C2)tJo vrc2 

Again it is seen how wq = is the physical position of the resonance, the point at which 
the scattering volume changes sign. For cjq < the system has a low-lying bound state 
and as wq changes sign, the bound state turns into a resonance. 

6.2 Properties of p- wave Feshbach resonances 

In the finite density atomic gas there will be three length scales, and thus it is 
possible to construct two dimensionless parameters. The interparticle separation is of 
order l/y^mep, a second length scale is provided by rn^g^ , while a third is the range of 
the forces Re ~ 1/A. The first dimensionless parameter is C2 defined in Eq. (6.11). The 
second is 

7p = vn^'^g^^. (6.17) 

In order for the results to be independent of the precise nature of short-distance physics, 
the interparticle separation must be greater than Re- This translates into the require- 
ment 

7p < C2. (6.18) 



The value of 7p determines the width of the Feshbach resonance in the same 
manner as 7s determines the width of s-wave resonances. The main difference between 
these is that while 7s increases with decreasing density, the opposite is true for 7p. Thus, 
l^wave resonances may be made arbitrarily narrow by reducing the particle density and 
it is possible to construct a pertubative theory, accurate across the Feshbach resonance, 
for densities achievable in current experiments. 

That the p-wave resonances are intrinsically narrow comes from the presence of 
the centrifugal barrier, which adds a long ranged 1/r^ tail to the interatomic potential, 
as in Ucs{r) = U{r) + l/m,j.r^. The width of a low lying resonance may then be 
estimated by computing the decay rate through the effective potential and is found to 
be dominated by the presence of the long-ranged tail [34] . 

Note how, even after the shift to a physical detuning, the scattering amplitude 
(6.15) still depends on the regularization scheme through the parameter C2. If C2 <C 1 
then the ultraviolet cut-off indeed drops out of the problem. However, if C2 ^ 1 then 
fco ~ —A and the bare molecule may be integrated out, resulting in a p-wave single 
channel model. Feshbach resonances with C2 ^ 1 are called weak while those with 
C2 ^ 1 are termed the strong resonances. Unlike 7p, C2 is an intrinsic property of the 
Feshbach resonance in question and depends only on the physics of the resonance, it 
cannot be changed simply by changing the density. 

Both of the parameters C2 and 7p control the perturbative expansion in powers of 
the coupling constant g. The usual approach is to apply mean field theory to the system 
described by the Hamiltonian (6.1). This approach has been used in Refs. [35, 17, 34]. 
The mean field approach relies on the smallness of the coupling g which means that both 
7p and C2 must be small. This is not necessarily the case since C2 is an uncontrollable 
parameter. 

The breakdown of mean field theory for strong p-wave resonances does not imply 
that no other techniques exist with which to investigate the system. Fluctuational 
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corrections arise from two distinct regimes of momenta, of the order of the inverse 
interparticle spacing and of the order of the ultraviolet cut-off. Corrections to the 
many-body problem are governed by the former and are small as long as 7p is small. 
At momenta of the order of the ultraviolet cut-off only virtual particles propagate and 
thus corrections to mean field theory from large momenta must be few-body, resulting 
in a renormalization of the few-body coupling constants. 

The situation is quite similar to the wide resonance s-wave superfluid studied 
in chapter 3. Here, the interactions between fermions were strong (and the scattering 
problems studied dominated by momenta of order a^^) and led to non-trivial renormal- 
izations of the few-body coupling constants, in particular it was seen how O)-, ~ 0.60. 
However, the many-body physics was dominated by physics at small momenta (p <^ a~^) 
in a perturbative expansion in the gas parameter. As long as the gas parameter was 
small, the mean field theory remained essentially valid, with the few-body couplings 
non-trivially renormalized. 

Finally, it should be mentioned that the p-wscve Feshbach resonance investigated 
in the experiment [30] is narrow and strong. This estimate was performed in Ref. [34]. 
The existing mean field theory is valid for the narrow and weak Feshbach resonances, 
and is in principle not valid for this particular Feshbach resonance. The few-body 
physics leads to a non-trivial renormalization of coupling constants, and indeed it also 
leads to bound states of three identical fermionic atoms interacting close to a strong 
Feshbach resonance [49, 41]. A reliable method with which to investigate the strongly 
resonant p-wave condensates is thus needed and this is the subject of the next chapter. 



Chapter 7 
Strongly resonant p-wave condensates 



While the weak and narrow Feshbach resonances are amenable to a mean field 
treatment [35, 17, 34], it is the strong and narrow p-wave Feshbach resonances which 
are most relevant to current experiments and properties of these are calculated in this 
chapter. The main results of this chapter were presented in Ref. [49] . 

An interesting effect in the vicinity of a strong resonance is the appearance of 
bound molecule- atom (trimer) states with angular momentum 1. This effect was first 
noticed by Y. Castin and collaborators [16]. These trimer states are quite unusual. 
They are very strongly bound, with a size of the order of the short-disance physics. 
Re, the same as the closed- channel molecule. The binding energy is very large (of the 
order of A^/m) and only weakly dependent on detuning from the resonance. Below, the 
binding energy of the trimer as a function of the strength of the resonance is found. In 
particular, the critical value of the strength of the resonance at which the trimer state 
appears is calculated. 

The two-channel Hamiltonian (6.1) is valid both for weak and strong resonances. 
In a manner similar to the s-wave case (see chapter 3) it reduces to a one-channel model 
in the limit of infinitely strong resonces, i.e. as C2 -^ cxd. However, it is interesting 
to investigate the manner in which the physics evolves under the strengthening of the 
Feshbach resonance, in particular the appearance of three-body bound states for strong 
resonances. Thus the starting point will be the two-channel Hamiltonian. 
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As discussed in the previous chapter, in order to study the two-channel Hamil- 
tonian (6.1), the couphng constant between the open and closed channels needs to re- 
flect the suppression of the molecule wave-function outside the range of short-distance 
physics. Different regularization schemes are possible. The two choices of cut-off func- 
tions used below give results which are in qualitative agreement, although they do not 
completely agree quantitatively. Arguments will be given as to why the results are still 
qualitatively trustworthy. 

A recent experiment [30] studied a gas of p-wave Feshbach molecules in '^^K. 
Unfortunately, the gas was found to be quite short-lived, with a lifetime of about 2 
ms. Ref. [30] studied whether the short lifetime might be due to dipolar relaxation but 
found that the lifetime was shorter than predicted from losses due to dipolar relaxation 
alone. The presence of the trimer state opens up for the possibility of inelastic collisions 
in which the scattering of two molecules results in an atom and a trimer with large 
kinetic energies. The lifetime of the p-wscve superfluid due to these decay processes is 
estimated below. 

One question which is not addressed below is the nature of the true ground state of 
the system. In order to properly answer this question, the molecule-molecule scattering 
problem needs to be studied and the few-body coupling constants computed, as in Ref. 
[34] for weak resonances. A diagrammatic technique similar to the one studied in section 
3.6 may be employed, but the solution of the resulting set of integral equations will be 
technically very difficult. 

In this chapter, first the diagrammatical expansion used to study few-body prob- 
lems will be discussed. Next, the three-body problem will be solved for different choices 
of the ultraviolet cut-off and it will be studied how the system evolves under the strength- 
ening of the Feshbach resonance. Finally, the stability of the gas will be considered. 
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1 ^ " 5 

Figure 7.1: The numbering of fermion propagators and vertices as explained in the text. 



7.1 The three-body problem 

In this section, the three-body problem consisting of three identical fermionic 
atoms interacting close to a Feshbach resonance is studied. It will be assumed that the 
system is tuned to the BEC side of the resonance, that is 

c^o = '-^ < 0, (7.1) 

1 + C2 

and thus a physical bound two-body state exists. 

It is convenient to explicitly write down the Feynman rules governing the dia- 
grammatic expansion. These may be derived from the Hamiltonian, Eq. (6.1). In the 
Hamiltonian, each vertex converting two atoms into a molecule (or vice versa) con- 
tains two fermionic operators. Thus the vertices commute and they may be arranged 
in any desired order. The atom-molecule scattering problem contains precisely one line 
of fermion propagators. -"^ Associate a number with each fermion line, from incoming 
to outgoing, and also a number with each vertex, as illustrated in Fig. 7.1. Then the 
Feynman rules are read off from the Hamiltonian: 

(i) At vertex i where a molecule of spin n is created associate a factor of {pi+i —pi)^/2 
where pi is the numbered fermion momentum. At vertex i where a molecule is 
annihilated, associate a factor — (pj+i — pi)^/2. 



^ This is also true for the two-boson irreducible vertex F appearing in the molecule-molecule scat- 
tering problem of section 3.6. The Feynman rules derived here apply equally well to molecule-molecule 
scattering. 
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(ii) Fermion contractions give the same signs as in the s-wave problem. In the three- 
body problem this is a factor ( — 1)"'+-'^ where n is the number of loops. 

(iii) Each vertex gives a combinatorial factor of 2 arising from indistinguishability of 
the fermions. 

(iv) Each fermion line gives a factor iG and each molecule a factor of iD^^. 

(v) Finally, each vertex joining fermions of momenta pi and pj+i gives rise to a factor 

igilPi-Pi+ilf^) = ig^iWi -Pi+il/2A). 

The atom-molecule scattering problem is studied in a similar manner to the cor- 
responding s-wave problem of section 3.5. The kinematics is chosen as follows: An 
incoming molecule of spin fi has four-momentum (0,E] and the incoming fermion 
(0, 0). The outgoing molecule has spin i^ and four-momentum {p,po + E) and the out- 
going atom (— p, — Po)- To calculate the s-wave scattering length, the total energy E 
reduces to a function ko('^o) ^ (this ensures that the molecular propagator is on-shell 
and will be defined precisely below) and to look for bound states the total energy must 
be taken less than kq. This is similar to the s-wave problem where E < Et^ = j. 

According to the Feynman rules above, the contribution to the scattering ampli- 
tude from the Born diagram of Fig. 7.2a has the value 

-iTl,]){p,po) = 2iZg%p,G{P,po + E)^ {\p\/K) ^ {\p\/2K) . (7.2) 

For proper normalization of the external propagators in the scattering amplitude, the 
residue Z of the molecular propagator at its pole is included (as in the s-wave three-body 
problem, see section 3.5). The Born diagram is proportional to p^py and goes to zero 
as p^pq ^0 (except in the case of scattering at zero energy with vanishing detuning). 
The first loop diagram, shown in Fig. 7.2b, takes the value 

-iTJSiP.P^) = 8Zg^ J ^q^q^G{q,qo + E)G{-q,-qo)Cm/A)Cm/^^) 
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(a) (b) 

Figure 7.2: The Born diagram (a) and the first loop diagram (b) contributing to atom- 
molecule scattering. The molecular spin has been made explicit while the kinematics is 
explained in the text. 



xD{q, qo + E)G{p + q,po + qo + E) 

x{p + g/2)a(p/2 + q),( (|p + g72|/A) e {\m + 91/A) . (7.3) 

Repeated spin indices indicate a summation over the corresponding spin. The propa- 
gator of molecules has been defined as D^y{q, go) = Sfj.uD(q, qo) where it has been used 
that the molecular propagator is diagonal in spin indices (independent of the particular 
choice of regularization) . 

Using the Feynman rules above, it is seen that the generalization to an integral 
equation is 

T^u{p,Po) = -2Zg^G{p,po + E)p^p,^{\p\/A)^{\p\/2A) 
7^;^r/^a(g, Qo)D{q, qo + E)G{-q, -qo) 
xG{p + q,Po + qo + E){p + g/2)„(p/2 + q),^ {\p + q/2\/^) C (Ip72 + g1/A) . (7.4) 

Exactly as in the s-wave problem, it is possible to integrate over the loop frequency with 
the result qo = —cf' /2m. Let po = —p^ /2m to achieve the same frequency dependence in 
both T-matrices of the equation. Define the "on-shell" scattering amplitude T^y{'p) = 
T^v{p, —p'^/2m). In order to proceed, it is advantageous to scale out the cut-off by 
measuring momenta in units of the cut-off, A. Additionally, measure the total energy 
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in units of IS? jm. Then the integral equation is 

-^^ [ q'dq I ^T,^{q)D{q, -q'/2 + E)G{p + q, V/2 - q'/^ + E) 

x{p+ g/2)a(p/2 + g).e (ip + m\)^ i\m + ^i) • (7.5) 

In writing this equation, for simpUcity the dimensions of the propagators G and D have 
been extracted and the remaining G and D are dimensionless. The integration J dVt^ 
denotes an integration over directions of q. 

The tensor T^y{p) must have the general form 

T,,{V) = T\p)5,. + T\p)P^ ^ Yl T\pKAp), (7.6) 

^ 1=1,2 

with a set of basis tensors defined as u]^y{j)) = 5^j_u and v?nu{p) = PfiPu/p^- It is then 
possible to derive two coupled integral equations for the coefficients T^ and T^. Multiply 
Eq. (7.5) by w'^uip) and sum over spin indices. On the left hand side a product of basis 
tensors appears: 

KMn-^M = \ = U,.^ (7.7) 

V / ji 

The matrix U is clearly invertible. Upon further multiplying the integral equation by 

the inverse matrix it becomes 

T^{p) = -2Zg^mG{p,-p^/2 + E)p^52,j 
2m^g^k p 2 



-7r2 / 



(■00 

\ q^dqD{q,-q^l2-,E)aj,{p,q,E)T\q), (7.8) 

JO 



which is a set of two coupled integral equations. The matrix a is dimensionless and is 
given by 

^G{p + q, -p^/2 - q^l2 + E)u\Mn\M 
x(p + g/2),(p/2 + q),i (|p + ql2\) i {\pl2 + ^"j) . (7.9) 
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In order to proceed, it is necessary to specify the precise form of the cut-off function 
^, as this affects both the integration kernel, through the matrix a, and the molecular 
propagator. 

7.1.1 The hard cut-ofF 

The coupling constant is now assumed to be constant up to the ultraviolet cut-off 
after which it drops to zero. That is, 

g(\p\)^gQ(A-\p\), (7.10) 

where G is the usual step function, and properties of the gas calculated with this as- 
sumption. 

With the hard cut-off, the molecular propagator was found in Eq. (6.12) to be 

5„ 



D{q,qo) 



-'fjLV 



(1 + C2)(^o - gV4m -UJQ + i{))+^^{q^/Am - qo - z0)3/2 tan 



-1 A 



(7.11) 



The physical on-shell molecule in the presence of the medium must have a pole as its 
four-momentum vanishes [38]. That is, the molecule is affected by the detuning loq and 
should be evaluated at a four-momentum [q, go + I'^oi^o)) such that D^^ I 0, kq{ujq) 
0. From Eq. (6.12) it is seen that k.q(ujq) satisfies the implicit relation 



ujo = Ko{uJo) + -—^ —{-Ko{uJo)f^'^tan ^ ^ ^. (7.12) 

I + C2 A ^^—mK{)[u>Q) 

For Wo ^ — the function kq reduces to the physical detuning. 

The residue of the molecular propagator at its pole is also found from the molec- 
ular propagator. It takes the value 



Z 



1 C2K0 ?,C2^/-mKQ _i A 
1 + C2 T-K-, T tan 



2 h? jm — kq 2A ^—rriKQ 

and reduces to Z « 1/(1 -|- C2) as loq -^ 0. 



(7.13) 



97 
7.1.1.1 The simplest hard cut-ofF 

The solution of the integral equation (7.8) requires the determination of the matrix 
a, Eq. (7.9). Note that this matrix does not depend on the strength of the resonance, 
C2 . It is instructive to first consider the approximation in which the momenta appearing 
in the integral equation are cut off at the value A. Strictly speaking this approach is 
not correct, as the momenta going through the vertices are not treated in a manner 
invariant under Galilean transformations. Instead it corresponds to the approximation 

H\P + Q/2\)H\p/2 + g])^a\m{m mSq. (7.9). 

Using this approximation, the matrix a may be evaluated to give 



an{p,q,E) 



ai2{p,q,E) 



4p2 
(p2 _ pg _^ q2 _ E){p^ +pq + q^ — E) p^ — pq + q^ — E 
8p^q p"^ + pq + q"^ — E' 

6/ + p2(5g2 _ i2E) + 3(g2 - E){q^ - 2E) 

2^2p2g,2 



(2p2 + g^ - 2E){p'^ -pq + q^ - E){p^ +pq + q^ - E) p^ -pq + q^ - E 

8p3q3 p2 _|_ pg _|_ g2 _ £; ' 

-3^-2p2 + 3g2 ^(q^ -Ef -p'^E ^ p'^-pq + q^-E 

a2i{p,q,E) = -^ + -^ log^-^ T^ E^' 

4j,^ ^p-iq p^ -^ pq -^ q^ — jij 



2E - 2p^ - q^ 
a22{p,q,E) = p-3 



2pq{2p'^ + 3g2 - 3^) 

„2 _ „„ ^ „2 _ RT 

.(7.14) 



+ (p\2p^ + 4q''-5E)+3iq^-Ef)lor^' ^^ + ^' ^ 



' p"^ + pq + q"^ — E 
Note that since the total energy E <0 there are no poles in the angular integration. 

The integral equation (7.8) may then be solved at any given value of the strength 
of the Feshbach resonance, C2, the total energy, E, and the detuning uq. Fig. 7.3 
shows the functions T^{p) and T'^{p) in the limit of very strong resonances with a small 
detuning. The s-wave scattering length is determined from 

«bf = ^THo) (7.15) 

evaluated at £^ = k,q{loq). For the specific choice of parameters 

Obf « 3A"\ at C2 ^ oo, wo = -10"^A^/m. (7.16) 
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0.4 0.6 

p/A 



Figure 7.3: The functions T^{p) and T'^{p) in units of ^^ evaluated using the simplest 
hard cut-off. 



The scattering length only weakly depends on the detuning as will be discussed below. 
Bound trimer states are found as solutions to the homogenous integral equation 
at -E < Ko- One bound state is found in the problem, at 



E^ PS -0.04AVm, 



at C2 -^ oo, ujQ = —10 A /m 



(7.17) 



The binding energy is seen to be proportional to A^, thus this three-body bound state 
is of size of the order of Re- The method used to find the binding energy is the following 
(also discussed in chapter 3): A bound state in the problem shows up as a pole of the 
scattering amplitude at a negative total energy, the binding energy E^. This in turn 
corresponds to a solution of the homogenous integral equation at that energy, thus to 
an eigenvalue equal to 1 of the integration kernel. Fig. 7.4 shows the behavior of the 
largest eigenvalues as functions of the total energy E. It is seen that an eigenvalue equal 
to 1 is obtained for a total energy of E\^ ~ — 0.04A^/m. 



7.1.1.2 



The invariant hard cut-ofF 



It is now interesting to see if the simplified approach above indeed captured the 
correct physics. Therefore, the matrix a of Eq. (7.9) is now evaluated in an invariant 
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Figure 7.4: The largest eigenvalues of the integration kernel in Eq. (7.8) as a function 
of total energy E in units of Is? /m. The eigenvalues are evaluated at C2 — > oo and 
cjQ = — lO^^A^/m. The largest eigenvalue is seen to go through 1 a.t E "^ — 0.04A^/?tt,, 
corresponding to the appearance of the bound trimer. 



way, using 
i {\P + m\) e (l#/2 + g1) = e (A - p2 _ g2/4 _ p-. -) e (A - g2 _ p2/4 _ ^. ^ (7_i8) 

It is seen that the maximum possible value of the momenta is now 2A and the integral 
equation is solved on the interval p G [0, 2] . The matrix a may again be obtained exactly. 
For p, g < A it reduces to the result (7.14) above. For larger p,q the limits of angular 
integration depend on the cut-off and, although the integration is simple to perform, 
the resulting expressions are not very enlightening and will not be written here. 

In the limit of very strong resonances (c2 — > oo) the integral equation (7.8) reduces 
to 



T\p) = 
6 '■°° 



-2Zg^mG{p, -p^ /2 + E)p'^52,j 



7r2 , 
^ -/O 



1 



E - 3g2/4 - cjo + (3^2/4 - ^)3/2 tan-1 



-aj,{p,q,E)T\q).{l.l9) 



It is of interest to study the behavior of the physics as the strength of the resonance 
is increased. Had C2 been zero, the only diagram which would have contributed is the 
Born diagram whose value was computed in Eq. (7.2) above. The Born contribution 
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Figure 7.5: (a) The atom-molecule scattering length in units of A^-"^ and (b) the binding 
energy in units of A^/m, both as functions of C2. The position of the appearance of the 
bound state is seen to be at C2 crit ^ 3.3. 



does not have a pole at a negative energy and thus the bound state is absent for weak 
resonances. It follows that the bound state found for strong resonances must develop 
at a certain critical value of C2- This is indeed the case as depicted in Fig. 7.5. The 
scattering length Obf vanishes for C2 = since the Born term is purely proportional to 
p^Pu, see Eq. (7.2). Obf becomes negative for small values of C2 and diverges to — oo at 
the critical value 

C2,crit ~ 3.3. (7.20) 

The scattering length then reemerges from +00 to go asymptotically towards the value 
of 

Obf ~ 1.9A~\ (7.21) 

as C2 approaches 00. The critical value of C2 ~ 3.3 is exactly the value at which the 
bound state appears. For large values of C2 the binding energy saturates at 



^3 PS -O.llA^/m. 



(7.22) 



The values obtained here in the large C2 limit are somewhat different from the ones 
obtained above using the simplified hard cut-off, but the results are still qualitatively 
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consistent. 

7.1.2 Softening the cut-ofF 

The above resuhs may be tested by applying a different cut-off function. In this 
section, the coupling constant 5'(|p|) = g^i\p\) is chosen to be 

ff(|pl) = <7exp(V/A'), (7.23) 

which is more smooth than the hard cut-off used above. 

The first order of the day is to compute the propagator of molecules using the 
new cut-off. To this end, the fermionic loop S is needed and according to Eq. (6.10) it 
becomes 

g^ f , g4exp(-2gVA2) 



V { \^^^ - '^mK{q^/Am - qo - iO) 



Stt^ J Po — q^ /m — p2/4 -|- zO 
mg'^ I \/27r 3 ^/2tt 



37r2 ^" I 16 4 



=V |-c'i + 4(^2/4^ - (70 - iO) - c'2^((7V4m - go - i^f'^F ( 



A^/tti 

(7.24) 



The coefficients c'^ and Cg are the analogues of ci and C2 calculated in chapter 6 for the 
hard cut-off. They are slightly changed and take the values 

, _ mff^A^ , _ mVA 

''i"l2(2vr)3/2' '^2-3^2vr)3/2- ^^'^^^ 

The function F replaces the arctan present in the hard cut-off propagator and is given 

by 

F{x) = \/27rexp(2x)Erfc\/2x. (7.26) 

Erfc is the complementary error function 

Erfc(x) = 1 - ^ / e-^^dt. (7.27) 

V^r Jo 
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(a) 




Figure 7.6: The functions T^{p) and T'^{p) in units of ^^ evaluated using the "soft'' 
cut-off. 



Note that F{x)x^q ~ \/27r while F{x)x^oo ~ x ^1"^ . 

Using Eq. (6.8) the molecular propagator becomes 



D^,u{q,qo) 



-'flU 



(1 + c',){qo - qy4m - u^o) + c'^^{qV^^ " 'Zo - i^?'^F ["^^f^ 



(7.28) 



As above, wq = -. , / . The on-shell molecular propagators must be evaluated at 



D{q,qo + kq {^o)) with kq satisfying the implicit relation 



Wo = Ko{uJo) + 



m 



1 + cL A 



-KO 



(.o))^/^f(-!^). 



(7.29) 



The residue of the molecular propagator at its pole takes the value 



Z 



, 2c2mKo(t^o) / /3 2mKo{uJo)\ I -mKo{uJo) „ / mKo(a;o) 



TT 



-1 



(7.30) 



and reduces to Z ~ 1/(1 + 02) ^^ ^^^ detuning vanishes. 

In the limit c'2 -^ 00 and loq ^ the scattering amplitudes T^{p) and T'^{p) are 
shown in Fig. 7.6. These are obtained by solving the integral equation (7.8) using the 
exponential cut-off. The functions are seen to approximately match those obtained in 
the simplified hard cut-off problem, see Fig. 7.3. 
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Figure 7.7: Behavior of the scattering length in units of A ^ under a change of strength 
of the resonance. 



Fig. 7.7 depicts the scattering length as a function of the strength of the res- 
onance. The detuning has here been taken to vanish. The critical strength of the 
resonance, above which the bound trimer state is present, is found to be 



C2,crit ~ 4.6. 



(7.31) 



The results presented above for the hard and exponential cut-offs do not match 
completely, yet the qualitative behavior is found to be the same. It should be emphasized 
that the binding energy for strong resonances goes as IS? jra and is very large. The 
binding energy is thus sensitive to the physics at large momenta, i.e. to the specific 
choice of the regulator 5'(|p|). The question is then, does the existence of the trimer 
depend on the choices of fi'dp]) presented above, that is, is it an artifact of the model 
and not the actual physical phenomena? To answer this question, consider the system 
close to the formation of the trimer, i.e. C2 close to C2,crit- Here, the binding energy of 
the trimer may be arbitrarily small and is a low energy phenomenon, insensitive to the 
choice of regulator (on the other hand, the precise value of C2,crit) being proportional 
to A, does depend on the choice of 5(|p|)). As C2 is increased from C2,crit) the binding 
energy E-j, quickly becomes large. However, once the trimer has formed, changing the 
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regularization scheme to make the model more realistic at short distances does not make 
the trimer disappear. Indeed, in order for the trimer to disappear at some c^ > C2,crit5 
the binding energy would have to first become large and then again drop to zero. The 
disappearance of the trimer would then again be a low energy phenomenon and would 
be predicted by the theory, contrary to what was found above. In this sense, even 
though the binding energy of the trimer is large, the existence of the trimer is a result 
of low energy physics. In conclusion, although the precise value of the binding energy 
of the trimer depends on the choice of regulator, the presence of the trimer state for 
strong resonances does not. 

7.2 Stability of the p-wave superfluid 

It is important to estimate the lifetime of the p-wave superfluid. In the strongly 
resonant BEC regime studied here, one of the main decay channels will be inelastic 
collisions in which two molecules turn into an atom and a trimer. In the process, the 
large binding energy of the trimer will be released as kinetic energy and the particles 
are lost from the system. 

As in the case of s-wave collisonal losses studied in section 5.2.1 the loss rate will 
be written as 

nb = -ainul. (7.32) 

rib is again the density of bosonic molecules. The inelastic relaxation rate is given by 

2 

Oil! = — {kiain{ki)) . (7.33) 

m 

Here, ki is defined as the relative momentum of the incident molecules and the average 
is taken over ki. The inelastic scattering cross section is given by [44] 

cJin= I \fin?d^^. (7.34) 

/in is the inelastic scattering amplitude into the final relative momentum kf and is of 
order i?e, just like the three-body scattering length above was found to be proportional 
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Figure 7.8: The processes which lead to the inelastic losses, Eq. (7.33). The trimer is 
illustrated as the binding of a fermion and a molecule. 



to 1/A ~ Rg. This may be seen from estimating the behavior of diagrams of the type 
shown in Fig. 7.8. These diagrams are similar to the form studied in the s-wave four- 
fermion problem but in the p-wave formalism. The final relative momentum arises from 
the kinetic energy released in the binding of three fermionic atoms into a trimer, thus 
/jj ~ A ~ 1/Re and 

(7.35) 



ttin ~ • 

m 

Additionally, there are inelastic decay channels similar to those studied in section 
5.2.1 in which a molecule decays into a deeply bound molecular state in the presence 
of an additional atom or molecule. Estimates similar to the above show that the decay 
rate due to these processes are of the same order as Om. Note that these are present 
both for weak and strong resonances, unlike the decay to the trimer state. 

Consider the experiment described in Ref. [30] in which it was found that the 
lifetime of diatomic molecules was quite short, of the order of 2 ms. It is interesting to 
estimate whether this could be explained by the losses due to the inelastic processes. 
The decay rate is Fin ~ a[^n and is density dependent. The atomic density in the 
experiment was n ~ 7 x lO^-'^^cm^^ and Re estimated by the size of the molecules to be 
.02 gq. Thus the decay rate is of the order of 10 Hz, corresponding to a lifetime of about 
100 ms. However, the prefactor to the estimate is uncertain and it is also unknown 
whether the losses observed in Ref. [30] would depend on the density. In conclusion, it 
is uncertain whether the observed lifetime of diatomic molecules in Ref. [30] is limited 
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by inelastic collisions. 

Shortly after Ref. [49], containing the work described in this chapter, became 
available on-line, the paper [41] by Jona-Lasinio et al appeared in which similar conclu- 
sions were obtained for the p-wave superfluid. 



Chapter 8 



Conclusions and outlook 



In this thesis, the s- and p-wave paired superfluids have been investigated. In the 
wide-resonance s-wave crossover, properties of the BEC regime were computed in the 
gas parameter expansion. Higher order corrections were found, and the mass imbalanced 
BEC was considered. It was studied how p-wave resonances naturaUy fah into categories 
of weak and strong depending on the precise details of the Feshbach resonance, and it was 
demonstrated how the strongly resonant superfluid is unstable towards the formation 
of strongly bound trimers with angular momentum 1. The rate of losses arising from 
collisional relaxation was estimated and it was seen how the predicted lifetime of the 
gas was slightly overestimated compared with the relevant experiment. 

A prominent element throughout the thesis has been the constant interplay be- 
tween few- and many-body physics. This was observed in the s-wave BEC regime, 
where few-body physics determined the molecule-molecule scattering length Ob in terms 
of the atom-atom scattering length a. Once 05 was known, the many-body physics was 
derived in terms of this few-body parameter. The decay rate of the superfluid due to 
collisional relaxation in both the s- and the p-wave superfluids was also demonstrated 
to depend on few-body physics. 

Possible future directions of study in the s-wave case include finite temperature 
extensions of the work. A study of the excitation spectrum at momenta of the order of 
the inverse scattering length, for which the underlying structure of the composite bosons 
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is probed, might prove fruitful. Finally, it would be interesting to compute higher 
order corrections to the chemical potential and ground state energy which explicitly 
demonstrate the underlying fermionic structure of the BEC of composite bosons. 

An exciting prospect is to restrict the p-wave superfluid to two spatial dimensions. 
Here, the BCS phase has been shown to be topological and will support vortices with 
non-Abelian excitations [62, 35], both necessary ingredients for a fault tolerant quantum 
computer [42] . It remains to be seen whether a stable, two-dimensional j!>-wave superfluid 
can be achieved. 
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